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Thermal Convection from a Horizontal Wite 


By Lawrence J. Crane, Glasgow, Scotland!) 


The steady flow considered here is generated by the diffusion of heat from 
_ along thin wire, placed in a horizontal plane, into a uniform fluid which extends 
indefinitely above the plane. Uniform profiles are found for velocity and tem- 
perature in the case of a gas whose coefficients of viscosity and thermal conduc- 
tivity vary directly as the absolute temperature. The solution given is valid 
when the vertical distance above the wire is very much greater than the dia- 
meter of the wire. The solution is a good approximation to flow in the atmo- 
_ sphere so long as the vertical distances are not so great as to render invalid the 
_ assumption that the fluid is uniform throughout. The solution is valid even 
_ when the density differences.in the flow are large. 


a ee 


-_™ 


Sol. , Lia Oe 


Equations of Motion 


centre of the wire, x-axis vertically above the wire and y-axis in the horizontal 
plane in the direction perpendicular to the wire. The equations of momentum 
and continuity for steady laminar flow are 


Ox Oy 


Se ee Pee bea Siow. fh 0 Be 8 a a in iees aie ad 


respectively (the quantities 9) and mj, which are used later, refer to the sur- 
rounding atmosphere); and g is the acceleration due to gravity. 


_ gradient term has been omitted because the pressure is, in this case, approxi- 
_ mately constant in the mixing region. The last term in equation (1 ) is due to 
_the buoyancy of the fluid. 
The diffusion of heat through the fluid is governed by the equation: 


dy dy \* dy}? 


1) Mathematics Department, The Royal College of Science and Technology. 
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saavon 


A system of rectangular Cartesian coordinates is taken with origin at the’ 


ti) Ou 0 Ou 
Qus tern =a (u 5) —¢(@- a) (1) 
and 
(0 u) < 0(@ v) ia. (2) 


‘respectively, where (u,v) are the components of velocity with respect to, the 
—(%, y)-axes; 9 and yw are the density and coefficient of viscosity of the fluid 


Equation (1) is the boundary layer equation of momentum. The pressure _ 


Cy 7 » (nT) _ 2 ( a (3). 


we 
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where T, C, and x are the absolute temperature, specific heat at constant 
pressure end coefficient of heat conduction of the fluid respectively. (7) and 
%9 Tefer to the surrounding atmosphere.) Equation (3) is the boundary layer 
form. of the equation of diffusion of heat. 
Equation (2) is automatically satisfied when a stream function y defined by 


oy _ Oy 


Cy 


is introduced. Equations (1) and (3) are simplified by introducing the von Mises: 
change of variable from (x, y) to (x, y). Equations (1) and (3) then become: 


Ou 0 Ou & nies 4) 
‘Ox ~~ (Ow (wow 5] U (1 a (4) 
and 
1G, flame oT 
noe = oy (0 Gy) 3) 
respectively. 


In the case of a gas, in which x oc w oc T and in which C » 1S a constant, 
equations (4) and (5) simplify to 


a 0 Ou a 
Dy, ~ Ho Go Ow (u ) + U (6) 


and 


Op _ % Qo 9 (u *) 
Ox Cpe i Op)? 


respectively. The variable y in equation (7) is defined by 


(Ren (Ea ape 


(This last result follows from the fact that the pressure is approximately 
constant over the boundary layer.) 


Solution of Equations (6) and (7) 
A solution of equations (6) and (7) of the type 
u = x? f(x) , 
4 (8) 
p= x" h(n) , 


where 7 = y/x", is now sought. The functions / and h satisfy the following 
boundary conditions: 


Vol. X, 1959 Thermal Convection from a Horizontal Wire 455 


On the vertical axis 
n=0, f(y) =A (ny=0, 


and on the boundaries of the flow 


H=tm, f= o=—0. 


A third condition on the motion is obtained by integrating equation (7) between 
the limits y = 0 and y = y,— the boundary of the flow. This gives 


Vb 

cH 

ral pdy=0. (9) 
0 


It follows from (9) that the total flux of heat J across a horizontal plane is a 
constant of the flow which is in fact equal to 


vb 
2 [ pay. 

i) 
When the expressions for g and yp given by (8) are substituted in (9) it is found 
that m + = 0 and 

"b 
l= 2 | h(n) dn . 
i) 


Two more conditions on the indices m, and # are obtained when the expres- 
sions for u and @ are substituted in (6). These relations, which follow from the 
requirement that the buoyancy and viscous terms in (6) shall be of the same 
- order of magnitude when + is large, are 


p—1=2p—-—2n=m-—}>. 
Similarly the terms in (7) are of the same order of magnitude for large ~ when 
—l=p—2n. 


These relations determine the values of m, n and p as — 3/5, 3/5 and 1/5 respec- 
tively. These values of m, n, and p entail that 4 varies between finite limits. 
_ Equations (6) and (7) now reduce to the ordinary differential equations 
a? 3 CROTPR h 
is rates = eee ¥ mre 1 
gt ri Ho 00 Gy FP) + 5 (10) 
and 


d ' 
— Sh tan FS (fh') , (11) 
Dp 
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respectively. Equation (11) can be integrated to give 


— 5 nh= "20 fh, (12) 


when the boundary condition h’(0) = 0 is used. 
The following changes of variables are now made 


a 2 n 
pt F(é), Peso He) cand = Se 
j= FO ag 0? HEE) iaae, 
where « is a certain constant. Equations (10) and (12) now become 
J60 
6F—18§F’=5- ge (FP) +16 = (13) 
and : 
—1866H=5FH, (14) 


respectively. (o is the Prandtl number of the atmosphere Cy f49/% .) The bound- 
ary conditions are: 
On the axis of the flow, i. e. on € = 0, 


EO se 0, 


and on the boundaries of the flow which may, by choosing a suitable value of «, 
be taken at € = + 1%), P=H=0 


For the special case of o = 5/9, equations (13) and (14) have the following 
exact solution which satisfies the boundary conditions . 
Peet 2. (15) 
Determination of an Approximate Solution of (12) and (13) 
when @ Is Close to 5/9 


Since the value of o for air is about 0-733, it is of interest to extend the 
solution given above to cover values of o near to 5/9. 
The method used to solve equations (12) and (13) consists in expanding F 


and # in the series: 
FoRt PRR +, 
H=H,+6H,+ PH,4+---, | 
where fy = Hy = 1— & and B = 9/5 —1. (In the case of air 6 is about 0-32.) 


(16) 


*) When the limits of & are fixed the value of x follows at once from the heat flux integral. It is 
in fact given by 
9 & ee |" 
2 \18 fo Qo 
when the boundaries of the flow are taken as E= - 1. 
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; The equations for F, and H, are then found by substituting the series (16) 
in (12) and (13) and equating the coefficients of 6. The equation for F; is 
Ti 12 9 Ub 18 2 16 ‘ 32 
Q-MA-Sea-wRa(S a-a+sla- Fe. 


This equation reduces to the following equation of hypergeometric type: 


=I 4 er Ts. 
z(l—aer+ {e+ z Lae ath pe 


5 10 ( (17) 


when the changes of variable d&/dz = €(z)/(1 — z) and z = & are made. The 


boundary conditions on ¢(z) are 
at2—=0,-¢ is finite andatz=1, C=O. 
The solution of equation (17) which satisfies the former condition is 
C= 2-H (a; bee) 2)— 1, 


where the usual notation for the hypergeometric function is used and where 


3 4 3 
Saae Gh and Oia 


a+b= 
The constant of integration A is fixed by the boundary condition at z= 1. It 
is in fact 
I (c — a) I (c — b) 
“E (6 = Ga DT (Cc) 


Finally 


the lower limit is unity because , vanishes on the boundaries of the flow. 
_ The function H, satisfies the equation 


Z&(H, + H,) +F A, +h H, =0. 


This equation can be integrated to give 


z 


/ aye hogar]. 


0 


H, = (1—2) 


The constant C is determined by the condition 


a 
[Hid=0. 
0 
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This condition follows from the fact that the total flux of heat from the wire 


is constant. 


Change of Variable 


The coordinate y is related to € by 


dy dé J on? 
yf B=o ol fe 35° ge H 


"dé ee a? 
=o elcrey + [teehee H]. 


0 


_ (18) 


When o = 5/9 (18) becomes 


y= 6 x70 lee [tanh 1Eé4 Tae as eae |] 3 


2 2 
= 6 x(Gr)-14 anh a a é]. 


Gr is the Grashof number of the flow; its value, for this type of motion, is 


125 ee ae)" ‘ 


Mo 


Conclusion 


A solution of the boundary-layer equations is found for the thermal flow 
due to a long thin horizontal wire which is heated. When the density of the fluid 
is sensibly constant, the motion is approximately of ‘similar’ type. In this case 
the velocity and temperature distributions are 


1 
a a x5 Fé) 
and 
Z a? 


oh 
pet eer 135° og x85 H(é) , 


where & is defined by 


“dg y % Oo y 
= we 1/4 
Ee 6 72/5 |/ Mo Bip (Gr)'4, 


When o = 5/9, these formulae simplify to 


om wild Fi y “ 0 cers 
Nic ae o6| we ef at at Y 
Zo sech. ls oon \ aa | =e sech ls : (Gr)] 
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1-2 


1-0 Caxis/Oo 


7. 
‘Similar’ 
0-8 - profile for 


Figure 1 
Plot of dimensionless velocity F against y/6 x2/5 Vx Qo/Mo for flow in air (o = 0-733) when 
Oaxis!Qo = 1/2, 2/38, 1. The ‘similar’ profile when o = 5/9 is also given. 


I-2 


profile for 


She g-2 
Bes 
e = 
% 0-4 08 12 16 2-0 
y, & Oo 
6x2) V Ho 
Figure 2 . 
Plot of e BBY fae ete 
73 we a To 
_ against y/6 x?!% Vu Qo! fo for flow 


in air (o = 0-733) when Qg»;5/Q9 = 1/2, 2/3, 1. The ‘similar’ pro- 
file when o = 5/9 is also given. 
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an Mi Og BX ae 
= 1 Tay eager Wem te ya, 

I a2 2 gy 1/4 
Ue ee g Pos sech E x (Gr) 


These formulae show that the width of the mixing region at a height % above 

the source is of order «(Gr)-1/4. This result agrees with the work of POHLHAU- 

sEN) for the width of the boundary-layer of the thermal flow on a vertical plate. 

When the density differences in the motion are appreciable the flow is no 

longer ‘similar’. This is because the second term in the square brackets cannot 
now be ignored. [The term referred to, namely 

z a 
ELY Seas 


H, 


is equal to (0,/o — 1).] As Figures 1 and 2 show, the greater the density dif- 
ferences the greater the width of the mixing region. For example the half-width 
of the mixing region when 0,,;;/Q) = 1/2 is about double the half-width of the 
mixing region when Q,,;,/09 is negligible. Now the density differences depend 
on the flux of heat convected from the wire. Thus the greater the flux of heat 
from the wire the greater the width of the mixing region, provided other 
conditions remain unchanged. 


0-4 0-5 0-6 0-7 0-8 0-9 1-0 


F(z) 0-610 | 0-563 | 0-514 | 0-463] 0-410] 0-354] 0-295] 0-232] 0-163] 0-087 0) 


A, (2) 0-094 | 0-055 | 0-008 |-0-025 |-0-052 |-0-071 |-0-083 |-0-082 |-0-071 |- 0-047) 0 


The data presented in Figures 1 and 2 are obtained by approximating to F 
and H by the first two terms of the series (16) respectively. When o = 0-733, the 
error is expected to be only a few per cent. This estimate is based on the assump- 
tion that /, and H, and subsequent terms are of the same order of magnitude 
as /, and H, on the axis of the flow. 


I am grateful to Professor D. C. Pack for the interest he has taken in this 
work. 


Zusammenfassung 


Eine exakte Lésung fiir die natiirliche Konvektion von einem langen horizon- 
talen Draht wird gegeben, selbst fiir den Fall grosser Unterschiede in der Dichte. 


Wenn die Unterschiede in der Dichte klein sind, dann ist die Strémung eine «ahn- 
liche». 


(Received: February 2, 1959.) 


8) E. PoHLHAUSEN in a paper by E. Scumipt and W, BrcKmANN, Forsch. Ing.-Wes. 1, 391 


(1930). — This paper is quoted by ScuiicuTiNG, Boundary Layer Theory, 1st ed. (Pergamon Press, 
London 1955). 
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A Quasi-Iterative Method 
for Computing Interface Temperature Distributions’) 


By FREDERICK F, Line, Troy, New York, USA2) 


1. Introduction 


In the study of friction and wear, the latter case especially, it is of impor- 
tance to know the interface temperature distribution of the surfaces in sliding 
contact. Such calculation is known to be difficult. In what follows, a simple 
model of bodies in sliding contact is proposed for the experimentalists and a 
relatively simple method is developed for the computation of the interface 
temperature. The method is applicable to more complex models. 

Given a stationary semi-infinite solid with a finite plane area (21x 25) 
protruding slightly above a plane surface of the solid. When another semi- 
infinite solid is slid steadily across this protruding area, temperature rise due 
to friction results. If radiation loss across the boundary is neglected, the 
problem of finding the interface temperature rise can be formulated as follows: 

The steady state temperature u(x, y, z) for a point in a semi-infinite solid 
(constant thermal conductivity k,) with zero imitial temperature, and zero 
temperature gradient normal to the boundary except over a finite area 
(-lsx </,-—b=y <b,z=0), where the rate of heat input per unit area 
is g,(x, y), can be expressed as(1}) 


1 6 
1 a dé d 
u(x, NK; z) = pee / / q(é, ”) et, (1) 
She Te 


where 
r= [le 8 + (y — 9) + ape. 


The quasi-stationary temperature v(x, y,—z) for a semi-infinite solid 


_ (constant thermal conductivity K, and diffusivity #,) with zero initial temper- 
ature, and zero temperature gradient normal to the boundary except over a 


- 


- 


» 


= 


finite area (-l Sx S<1,-—b Sy S),z2=0), which moves with a uniform 


1) The research was supported by the United States Air Force Office of Scientific Research, 


Air Research and Development Command, under Contract No. AF 49(638)-67. 


2) Rensselaer Polytechnic Institute, Department of Mechanics. 
3) Numbers in brackets refer to References, page 473. 
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velocity V in the direction of x, can be expressed as [1] 


a4 sd he iT ds dy (2) 
v(x, ¥, 2) = nie i [ale m) € ae 


—l —b 


where x here denotes the Lagrangian coordinate with respect to the centroid of 
the moving heat source, ,(x, y) is the rate of heat input per unit area. 
If the heat generated at the interface is g,(x, y) then 


q(x, ¥) + dalx, ¥) = Go(*, ¥) (3) 


and the interface temperature is usually expressed by the equivalence of 
equations (1) and (2) for -1 <x S/,-b Sy Sbandz=0 


u(x; 9, 0) = v(%,,7,.0) x (4) 


From equations (3) and (4),'a singular integral equation in the unknown 
qi(x, y) is resulted. Once the inversion is obtained, i. e. 9,(x, y) found, either 
equation (1) or (2) can be used to compute the interface temperature. 

Numerical methods have'been applied''to similar problems. Also trial and 
error schemes have been used. to find q(x, y). By far most of the approximation 
methods employed have sprung from BLox’s original work [2] which utilizes 
the following idea: For qg, independent of x:and y, assume q, and q, are inde- 
pendent of x and y. Integrate equations (1) and'(2) to give u(x, 0,0) and 
v(x, 0, 0) in terms of 9g, (¢2 = Uo — 1). Then arbitrarily equating u(x, 0, 0) nar = 
v(x, 0, O) max; this provides an additional equation for the unknown q,. Using 
the value of g, thus found, w*(x, y, 0) and v*(x, y, 0) are computed, where u* 
and v* refer to u and v for g, and g, equal to constants. The probable interface 
temperature, BLOK argued, must lie between w* and v*, and therefore can be 
sketched in. This leaves much to be desired, especially if a more accurate 
estimate of the interface temperature distribution is sought. 

In section 2 A an iterative method is outlined for the computation of the 
interface temperature. Section 2B gives a quasi-iterative method which 
furnishes the advantage of a faster convergence of the process than that 
outlined in section 2 A. Section 2 C shows the exact solution to the problem 
for the asymptotic case R < 1, where R is the Peclet number associated with 
the problem. Moreover, the iterative process of section 2 A is further justified 
in a quantitative manner for R <1. Section 2 D shows the asymptotic case 
R > 1, where the iterative process takes on a simpler form than the general 
case outlined in 2 A. Section 2 E describes a sample calculation. 

Section 3 shows that the interface temperature problem can be used as an 


analogue problem for inverting a class of singular integral equations of the 
first kind. 
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2. Method of Solution for |< b 
A. Iterative Method 


In subjecting similar problems to digital computation, CHao and TriccER{3] 
have found that for / < 6 and q independent of x and y the interface temper- 
ature is essentially independent of y (the exception being near y = + b). This: 
suggests that g, and q, are essentially independent of y for this case. In what 
follows it is shown that, based on the assumption that g, and q, are independent - 
of y, u(x, 0,0) = v(x, 0,0) can be evaluated analytically for 1<b and q 
independent of x and y. For simplicity, let K, = K, = K. Also let g, = gy m(x). 
From equation (3), ¢, = dg{1 — m(x)]. Define f=a K u(x, 0, 0)/qol, s = x/l, 
t= E/l, € =y/l and B = djl. 

Equation (1) becomes 


f(s) = és mi(t) F(s, t) dt, (5) 


where 


F(s, t) =InB law [a (254) | -nl ef, 
and 
F(s,t)~In2B—I\In|s—#| for 1, 


Furthermore, define g = 2 K v(x, 0, 0)/q9/ and R=v1/2a, equations (2) 


"and (3) yield 


1 
gls) = [ ult) Gls.) dt, ple) = 1 — mle), (6) 
—1 
where 
G(s, t) = e®&-9 T (s — t;'R).. 
The function 
P RVG ge 
Vesta C 


I(s—t; R)= 


_ for certain values of the parameter R have been tabulated in the literature [8]; 
and for large values of R, it can be approximated by.Hankel function of the 


owe Eh te a be 2 Ea ee es 6 bee Fes ee Bee 
: i 


zeroth order. oy 
Equating equations (5) and (6) as A ala at the interface for continuity 


[rend i= me () (Fl, + Gls, )) dt. (7) 
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Equation (7) is a singular integral equation of the first kind with a loga- 
rithmic singularity. The method of inverting this type of equation as outlined 
by CaRLEMAN [4] and Pearson [5] is impractical if one is interested in the 
solution for a wide range of the parameter R. For instance, PEARSON’s method 
calls for the expansion of the non-singular parts of the kernel in polynomials of 
(s — #); and he concluded that the time required to write out solution is about 
2 n3/2h, where vn is the order of higher of the polynomials representing the 
non-singular parts of the kernel. However, the existing method will be useful 
in inverting equation (5); the result can then be inserted in equation (6), which 
becomes a singular integral equation of the second kind. Thus f(s) = g(s) at 
the interface can be found by iteration. 

For B > 1, as stated above, F(s, t) = In2 #8 — In| s —#|. Equation (5) can 
be inverted [4] 


m(e) = —_ 2, [A Se 5 | Sane 


where 


1 ae Es 1 
_ =1 f Yi #zZ@) a 0 Lds 
L(z) : pas Exe) ——_ 
wt =f 

‘The problem of /(s) = g(s) for —1 < s <1 [equation (7)] can be resolved by 
the following process: 

Arbitrarily let mo(s), the initial trial value of m(s), be equal to zero, then 
the initial trial value of g(s) 


1 
gols) = | G(s, 4) at 
il 
and the resulting value of m(s) 


m,(s) = M(go) , 


where M(g,) is the linear integral operator [e i i 
abbas. quation (8)] operating o 
The first iterated value of g(s) (8)]] op & On £o(s) 


ex(s) = [ [1 — Mg] Gls, #) dt = eu(s) — THs), 


where 


: 1 
TO = f M(g,) G(s, #) dt 
rs | 
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after the m-th iteration 


where 


TH (go) es / M [To—) (2) ] G(s, t) dt . 


=! 


It may be seen from the kernels F(s, ¢) and G(s, ¢) that, although the former is 
symmetric and the latter skew-symmetric, the integral operator 


1 


P{z(s)] = / 2(t) F(s, t) dt 


is always greater than the operator 


_ for all values of —1 <s <1. The extent of the difference between P and Q 
depends on s and the parameter R. Of course P = Q when R= 0. 

In view of the above, g, = Q(1) and | M(g,) | = | P-(g,) |, P-1 being the 
inverse of P, is less than 1/(V1— s?) for —1 <s <1. The factor (1 — s)-12 
is associated with the root singularity at s = +1. In other words it may be 
demonstrated that | T"(g,)| < |g | for -—1 <s <1 and R + 0. Moreover, it 

_ may be demonstrated that 


| FO glib vin | T (gq) |< | go fe 


Thus equation (9) constitutes an alternating Neumann series whose convergence 
_ is assured. 


B. Quast-Iterative Method 


The rate of convergence of the above may be slow, in which case the 
following quasi-iterative scheme will always speed up the iterative process 
described in section 2 A such that no more than one quasi-iteration would be 
necessary in most cases. 

It is clear, that the closer the initial value of the trial function g*, the more 
- rapid would be the convergence process. Also, if there is an a priori method of 
finding a close-approximation of the initial trial function gf, the same method 


Vas 


— 
= 
-f 


— iP TY CU 
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can be used between every successive iterations. After BLoxk [2], define 


fins) = | Fis, t) dt [equation (5) with m = 1] 


and 
1 
g,(s) = / G(s, 4) dt [equation (6) with w= 1). 
= 


Then define 


— [&m(S) Imax 
Ler (ees 


By assuming m and ym to be constants and matching only the maximum 
temperatures of the mating bodies, m and « may be found. From equations (5) 
and (6) 

MU fin(S) Imax = (1 — ) [8m(S) Imax 


or 


m=—— and = 


Thus 


fol) = 4 tals) and enls) ==> Sel) 


must envelope, BLoK suggested, the true value of f(s) = g(s) for -1<s <1 
and the latter can be interpreted. 

Define gf, the initial trial value of g for the quasi-iterative method, to be 
(fp + g,)/2 arbitrarily knowing that this value would be much closer to the 
true value than g, in the last section. Using equation (8) 


my = M(go) 
and using equation (6) 
B*(s) = gy — TMs) , 
where, as before, 


1 
T" =/™ [T°] G(s, t) dt. 
eT 


Similar to equation (9) 


PIERS 2 Sig oO Sig Crab CoH Nh (10) 


% 
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Very approximately, T™ > A gy, T®) = A2 gy--- 7 = 4" g,. Therefore 


Note 


converges to 1/(1 + A). 
Define 


This first quasi-iterative value of g** just defined may be shown to be much 
closer to the true value than g*. Similarly a second quasi-iteration may be used 
to give an a posteriori check on the rate of convergence of the process. Of course 


C. Asymptotic Solution R <1 


For R <1, the problem of finding f(s) = g(s) [equation (7) ] is reducible to 


F(s) = | u(t [1+k(s—f)] (a—2kB—In|s—#)) de, (11) 
_ where 
F(s) = z { @—I|s—4) dt, «=In2, and —s 


1 


In the manner of PEARSON [5] equation (11) may be inverted exactly and 
~ equation (5) or (6) used to compute f(s) = g(s) for —1 Ss S1. To the first 


Tak 


etit? 


PR TA aN RPA 


BPS) a 


order of k, 


+ L(F’) 


1 A*(F) ale) 
1 +a 1+a 


| (Gaz -2) 
a 2 In2 We 
pa [Prete felt ase) — pg BA) 
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ue 
In2 p ne ya ba f A) pee Eee) 
beg ine ee ee i (eer es 


and the other integral operators are as defined in section 2 A. 

The complexity of equation (12), which is only accurate to the first order 
of k, justify the use of an approximate method in section 2 B. The fact that, 
for R <1, both equations (5) and (6) possess exact inverses makes it possible 
to study closer the convergence property of the iterative process outlined in 
section 2 A and 2 B. Equation (5) has an inverse m(s) = M(f), equation (8). It 
may be shown that equation (6) possesses an inverse for KR < 1 


p(s) = M(g) + N(g) , (13) 
where 
R Ppyie ene fer 
N(g) = = yi—s | Pon al + a)? (3) = te (g) 
6b — a ’ ’ 
— L(g) + 20528 peigy 4 74 Tene | 
Let . 


So= | G(s, t) dt=g+4+ &, 
at 


where g is the exact value of g(s) = f(s) and é, is the error in the initial trail. 
Using equation (8), 

MS aU a M (é») ? 
where m is the correct value of m associated with g. Using equation (6) 

J 
f=S+e = [1 —_m — M(e,)] G(s, #) dé 
PS) wt (14) 
=g— | Mle) G(s, dt, 
=) 


where é, is the error after the first iteration. 
Using equations (13) and (14) 


1 


c= — | (Mle) + Nleq) — N(eq)) G(s, ® at : 


—1 


= =e) + TM) , 
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where 
zk 


T(¢5) = / N(e) G(s, t) de. 
=] 
It is seen that N(g) is of the order of R and that it is demonstratably 
‘positive for —1 <s <1. Similarly 
eh Sede Te) 


Thus the convergence of the process is assured. Similar treatment can be made 
for the quasi-iterative case. 


D. Asymptotic Solution R>1 
The differential equation corresponding to equation (6) is 


02v 07u 02v .: Ou 
xt * Oy? + O28 — a” O(x/V) 


(15) 


recalling g corresponds to v(x, 0, 0). 
For R large BLox has found that 0?v/0x? and 0?v/0y? are essentially zero. 
Thus equation (15) is approximately 
0? 1 ov 


axe —mSC« Ox /'V) 6) 


The integral representation of the solution to the boundary-value problem 
governed by equation (16), which holds in a semi-infinite solid with zero 
initial temperature and zero temperature gradient normal to the pougesny 
except over a strip —/ < x S/, is 


2(s) =f od die naa ea eee ye (17) 


where the dimensionless quantities g, R etc. are as defined in section 2 A. 
Equation (17) is ABEL’s integral equation whose inverse is 


-1f a = (18) 
4 where k = 2 R/z?. 


-ZAMP X/30 
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For all practical cases, where equation (17) applies, k <1. Equation (5; 
may be written as 


jis) = [1 = ule) (@—In|s—#)) dt. (19 
Combining equations (18) and (19) for -l1 Ss S a. 


f(s) = fol8) — i nm “| («—In|s—¢)) de, (20 


where 


Equation (20) is a singular integral equation of the second kind; if k <1 anc 
the integral in equation (20) is bounded, the solution is expressible in the form 
of an alternating Neumann series. In other words, 


nw 


ta(s) = fols) + BD) (Br T(fo) , 


r=1 
where 
1 (Spry yar 
TH =| pe (n —In|s —#]) de. 
: f yi-t 
-1 -1 


E. Example 


An example is carried out for R = 1 and f = 4. The case of R = 1 is chose 
for illustration because on the one hand, the solution cannot be reached accu 
rately by the asymptotic solutions for R < 1 or R > 1, and on the other hanc 
the Kernel G(s, #) cannot be approximated by the Hankel function. Thus the 
case of k = 1 represents a case in the most difficult range of the parameter R 
The quasi-iterative procedure of section 2B is followed. The integrals ar 
evaluated with the aid of an IBM 650 digital computer. Integrations near th 
neighborhood of singular points are evaluated analytically while integratior 
over the non-singular portions are evaluated according to Srmpson’s rule. Th 
over-all computational accuracy is estimated at 11/,%/ at the worst locations 

Figure 1 shows gi = (fz + gg)/2, where f, and @p are as defined in sectio1 
2B. gf* and g*** are the result of the first and second quasi-iterations respec 


tively. Figure 2 shows a plot of fp, gz, and g***, It may be noted that f, and g, 
do not envelop g¥** at all points. 
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Figure 1 Figure 2 
Distribution of interface temperature for Temperature distribution g§** shown together 
R=1, B=4. — — g¥ first trial; with BLox’s bounds g® and /4. 
°0° g$* second trial; —— g*** final value. aS ye SS Fg 


3. Other Applications 


The interface temperature problem may be used as an analogue problem for 
_ inverting a class of singular integral equations of the first kind as a result of 
the above analysis. Given 


h(s) = / m(t) H(s, t) dt, (21) 


where h(s) is continuous and H(s, #) has an inherent singularity at s = ¢, a and 6 
are finite. The inversion of this type of equation have been discussed (see for 
instance [7]). The following method will be quite generally useful and simple 
to apply if (1) the kernel H(s, ¢) can be split into two parts, F(s, ¢) and G(s, 2), 
such that F(s, ¢) possesses an exact inverse; (2) numerical evaluation of a set 
_ of integrals are permitted; and (3) the particular H(s, ) are so split into F(s, ¢) 
and G(s, ¢) that certain convergence property is satisfied. 
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H(s, t) = F(s, t) + G(s, 2), = /uo F els) = | m(t) Gls, 4) dt, 


p(t) = P(h), the inverse of the equation 


b 
n(s) = | (0) F(s, t) dt 
Then equation (21) is equivalent to the following problem: 
f(s) = g(s), p(s) = mls) + ms) @Ss S08). (22) 


This equivalence is easily seen because equations (22) mean 


= [0 — m(t)| F(s, t) d = /m 


or 
b 


b 
IEO F(s, t) dt = [mo [Fis fe Gls Aig 


a 


and by the definition of p(t) and H(s, #), 


b 
= [ mle) H(s, # at 


To invert equation (22), let 


= | (0) — nO) G(s, 4 ae 


b b 


=| bl) G(s, 1) dt — [ ule) Gls, 0) at. 


But p(t) = P(h), w= P(f) and P(f) = P(g) as required by the problem, therefore 


— | Pie) G(s, at, (23) 
| | 
= | Pit) Gs,t 


where 
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Equation (23) can be solved by iteration. Of course the convergence property 
has to be examined in each case. Once g is found to any degree of accuracy 
desired 


Hs) = P(g); m(s) = p(s) — w(s) . 


Problems which could easily be treated in the above fashion, among others, 
are the Rijke-tube problem in fluid mechanics [6] and the welding problem in 
which the heat distribution of a travelling heat source is sought for a given 
temperature profile. 


4. Conclusion 


For the proposed model of sliding contact, temperature distribution of the 
interface due to friction may be calculated as accurately as desired. The quasi- 
iterative method reduces the number of necessary quasi-iteration to no more 
than one in most cases. An example is given to show this point. 

The method may be applied to more complex models and extensions can 
easily be made for the proposed model to accomodate a general heat input 
distribution due to friction. Moreover, the factor (1 — x?)/2 always occurs in the 
distribution function of the singular integral equation with logarithmic singu- 
larity; by a suitable transformation this factor may be cancelled [6]. This 
would simplify computation for those who prefer using digital computers. In 
other words the uncertainty at an end point experienced by CHAo and TRIGGER 
can be removed. 

Although the nature of the method is iterative, the labor involved is less than 
that required for numerical solution, given a degree of computational accuracy. 

By analogy the method may be used to invert a class of singular integral 
equations. 
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Zusammenfassung 


Es wird ein einfaches Modell fiir gleitenden Kontakt vorgeschlagen, zur Ermitt- 
lung der Verteilung des Verschleisses, die ihrerseits mit der Verteilung der Tempe- 
ratur lings der Gleitflache zusammenhangt. Dabei wird eine relativ einfache 
Methode entwickelt fiir die Berechnung der Verteilung der Gleitflachentemperatur. 
Die Methode ist im wesentlichen iterativ; im besondern verringert die hier ent- 
wickelte quasi-iterative Methode die Anzahl der Quasi-Iterationen betrachtlich — 
in den meisten Fallen auf nicht mehr als eine. Dieses Verhalten wird an einem 
Beispiel erlautert. Obwohl die dargestellte Methode iterativ ist, erfordert sie weniger 
Rechenarbeit, als fiir die numerische Lésung des Problems bei vorgeschriebener 
Genauigkeit notwendig ware. 

Es wird auch gezeigt, dass das Modell als Analogon fiir die Inversion einer 
Klasse von singularen Integralgleichungen erster Art dienen kann. 
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Energy Levels for the Coulomb Potential with Cut-off 


By Asp EL SADEK EL MeEticy and Manmoup AHMED EL SHERBINI, 
Alexandria, Egypt?) 


. 


Introduction 


The absorption band spectrum of alkali-halide crystals has been attributed 
to the excitation of electrons trapped in holes formed by the removal of negative 
ions. For the potential field of such electrons Trpss [3]?) proposed a Coulomb 
field with a horizontal cut-off which is represented by the potential function 


( 
29 


e2 


V= i - 
se ai ar (? eras i > Lor (Cte (1) 


where 7 is the radius of the hole and K, is the dielectric constant. He computed 
the electronic energy levels by numerical integration methods. WANNIER Ay 
however, has shown that the problem could be solved analytically, though not 
rigorously, for energy levels lying very high above the potential energy in the 
hole. The object of this work is to present a rigorous solution valid for low as 


well as high energy levels. This has been achieved by using exact expansions 
for the wave functions. 

1) Faculty of Science, University of Alexandria. 

) Numbers in brackets refer to References, page 477. 
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The Wave Functions 
The radial wave equation for the potential field (1) is 


d2F a? Ae ~ $ (t+ 1)\ 2. 
e ae aS F=0 (r<n), (2) 


@2E a 2a Trae = 
HHA EIEN ten, 6 


where « = m e?/K, h? and the energy E = — h? 2/2 m k?; k is our unknown 


which has to be determined. 
The solution of equation (2), which vanishes at the origin, is 


fms hose [fe (t~ ade)" “ 


and the solution of equation (3), which vanishes at infinity, can be written as 


=x fre ja(%) — x tan (z hk) 8x, 1+1/2(*) , (5) 

where 

x=2(2ar)"?, 
_ ot 
far y2(%) = 2% Ps ve (az) 7) 
and 
ba RAD = Oey #P oy 

Span) = — 29" aR pay erie (eR) ” 


% being the value of x at r = r, and the functions J, ,,(z) and N, ,,(z) being the 
confluent hypergeometric functions [1]. The functions / and F, represent the 
_ appropriate wave functions of the electron in the potential field (1). 

To evaluate the functions f, ,,(x) and g, ,,(x) as well as their derivatives it 
is necessary to derive their recurrence formulae. Applying the recurrence 
formulae given by PowELt [2] for the Coulomb wave functions, one can deduce 
the following formulae for the function f, ,,(x): 


(= 2) [lr 2) A] ane) = (m2 3) Mh ° 
+ 2mm (1—- =" =) f, ale), 
Pe fant) = [72 z ine ay] fam (9) (10) 


x (2m + 1)? 
2(2m + 1) fee A R2 | Fever - 


Similar formulae hold for the other function g, ,,(x) and its derivative. 
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Expansions of the Wave Functions 


For high energy levels, WANNIER has applied the zero energy limiting 
formulae, i.e. the asymptotic values when k tends to infinity, which are as 


Pa Fel) Toa A) (11) 
IP os —m+k 
ga) Re (2 Val)» (12) 


r(5 +m +k) 


where J2,, and Y,, are the Bessel functions of the first and second kind respec- 
tively. 

Now for exact representation, one can extend the result (11) to finite values 
of k by assuming the following expansion: 


fe, m(*) 7 Dey a, x” Jom+y(*) 3 (13) 
v=0 
On equating the coefficients of powers of « on both sides of this equation, it is 


found that a) = 1, a, = 0, a, = (2m + 1)/32 k?, ag = — 1/96 Fk? and in general, 
fory = 3, 


32 2 vya,=2(2m+yv—1)4,,—4,.,- (14) 


On the other hand, it has not been possible to extend the result (12) to 
finite values of k. However, one can express gp, ,,(%) in terms of f, ,,(x) functions 
by applying the expansion of N,_,,(z) given by MEticy [1]. Thus the exact 
representation of g, ,,(x) takes the following form: 


bsnl) = —hanl®) [5 v (a +m +8) + Fy (F—m + 2) | 


x2 
— (2m) —y (2m +1) + logf--| 


4% \2m x \Y 1 1 
5 oo wae >z>tajvth 
4 Gy es, v! (2) v)? : ae a 
a (15) 
akan Beare! = P(t ayo) GP 
r(5 +m+r) mas TI 
oe) (se) 1p? La ante 
Pa v! TS + 4 Tryya(X) » 
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where yp is the logarithmic derivative of the gamma function and the term 
v = 2 min the last series is to be discarded. The function x, m(X) could be easily 
calculated from (15) if one computes.the functions fr. yyo(2) from (13) and the 
recurrence formula (9). 


Evaluation of the Energy Levels 


The values of k& are determined by equating, at « = Xo, the logarithmic 
derivatives with respect to x of the functions and F, in equations (4) and 
(5). Thus, for s-states we have 

A fe rjo(*) Tr (1 iat, fr, 3/2(*) 


tanz k = at ee 
A’ 8, yo(%) + (1 — 2-4) & 5/0(#) (16) 


and for other states 
Priya) + B fy 1. 1/2(%) 
taitte ea Oe ae é 7 
8x,1-1)2") + B 8, 7. 1/2(*) a7) 


where 


8 81 —1/2\¥ 
= ycoty, B=1—-> Sian) 
x a Jis1)2%) 


x 1 eee We 


%» has been replaced by x for simplicity. Equations (16) and (17) do not contain 
the unknown & explicitly. However, their right-hand sides are slowly varying 
functions of k and so the equations can be solved for k by successive appro- 
ximations. Starting with k = » on the right, where x is the principal quantum 
number, two operations are sufficient to give an accurate value for k. For 
example, in the case where x = 2-62, the value of k for the 1s state is 1-161 
which is considered to be more accurate than the value given by Tipzs (1-155) 
or that given by WANNIER (1-153). 
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Zusammenfassung 


Die Absorptionsbanden der Alkalihalogenidkristalle werden auf Elektronen 


zuriickgefiihrt, die an leeren Gitterstellen eingefangen sind. Als Feld wird ein bei 
__kleinen Distanzen horizontal abgeschnittenes Coulomb-Feld angesetzt. Die Eigen- 


werte werden im ganzen Energiebereich aus einer genauen Reihenentwicklung der 


_ Lésungen der Schrédinger-Gleichung bestimmt. 


to i | 
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Analysis of Turbulent Boundary Layers with 
Zero Pressure Gradient 
By Joun R. Purrip, Pasadena, California, Wiss) 


1. Boundary; Layer Equations 


We are concerned here with the boundary layer developed along a flat plate 
of effectively zero thickness, with the leading edge of the plate at x = 0 and 
the plate lying along the x-axis. The free stream velocity of the incompressible 
fluid, U, is directed parallel to the x-axis. We denote by y the ordinate normal 
to the plate, and confine our attention to the region y = 0. We here assume 
that the boundary layer is turbulent right from the leading edge. 

Let w, v be the components of mean fluid velocity in the x and y-directions. 
If then, we equate for any volunie element of the fluid, dx-dy, the mean rate 
of gain of momentum in the x-direction to the x-resultant of the shear forces, 
and use the continuity requirement, 


Ou Ov 
ox + oy — 9° : eae 
we obtain 
Ou Ou i OT 
Pes Oyo. oy’ ee 


where Tt is the shear stress and @ is the fluid density. Equation (1.2) also follows 
immediately by analogy with the equations for the laminar boundary layer 


due to PRANDTI [1]*). Now, recognizing the turbulent and viscous components 
of the shear stress, we may write 


is a Ou 
ey => (kK ae v) Oy 5 (1.3) | 


where K is the eddy viscosity and y the kinematic viscosity. Now 


Ou 


2ae 


: (1.4) 


; 1) Address for 1958: California Institute of Technology. Permanent address: Commonwealth 
Scientific and Industrial Research Organization, Canberra, Australia. 
2) Supported in part by a grant from the Rockefeller Foundation RF 56160. 
3) Numbers in brackets refer to References, page 501. 
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where / is the mixing length as conceived by PRanptt [2]. In the present paper 
we introduce, and explore the consequences of, the hypothesis that 


l= ky", (1.5) 


where k and m are to be discussed in detail later. At this point we need only 
indicate that & is positive, non-zero, and a function of relative distance from the 
plate across the boundary layer and that m is a constant. 

(This introduction of a mixing-length formulation is quite empirical, and 
finds justification solely in the unified treatment, and the agreement with 
experiment, we are able to present later in this paper, and in subsequent work. 
Note that we could have introduced as a general expression for the turbulent 
shear stress 
B2 2m | Ou |n Ou 

“  [ay} oy ° 


_ The experimental fact that local drag coefficients over rough surfaces are 


independent of U leads back to exactly the case we consider here, 1. e., m = 1.) 
We now neglect v as small compared to K. This evades the problem of the 
laminar sublayer, and, also, as we see later, is inadmissible at the extreme outer 


edge of the boundary layer. Then use of equations (1.3), (1.4) and (1.5) in 


equation (1.2) leads to 


(1.6) 


- Ou/0y = 0 in the region we consider, so that we effect a typographical economy 
_ by writing 0u/dy for | 0u/dy | in equation (1.6) and subsequent equations. Here, 


and subsequently, we could reduce the number of equations by one by intro- 


- ducing a stream function. However, the treatment would become more compli- 


ors 


ee ee RO Nee Bt 


cated in other ways, so we retain equations (1.1) and (1.6) as the boundary layer 
equations. 
2. Reduction of Boundary Layer Equations to Ordinary Form 


Substitutions (2.1) used in equations (1.6) and (1.1) enable their reduction 
to the ordinary equations (2.2) and (2.3): 


oe ne fay xO) vy Un-2, Cav 2-282) Y-1 | (2.1) 
Si ielas cas =A dy dn uf) age (2.2) 
dl ite ee. (2.3) 
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This reduction is permissible only if the boundary conditions governing 
equations (1.6) and (1.1) can also be reduced by substitutions (2.1) into forms 
independent of x and y. This would be so for the boundary conditions of the 
corresponding laminar problem, if these could be carried over legitimately into 
this turbulent problem. These conditions, which we shall not use as they stand, 
are: 


Lt te A ey eee (2.4) 
KO = 6, (259) 
yoo, u>U, Le ae (2.6) 


Oy 


Use of conditions (2.4) and (2.6) would imply that equation (1.6) holds 
throughout the infinite region y = 0. In fact, there are two regions in which 
equation (1.6) fails. In this section we discuss the region close to the plate, 
deferring the question of the outer limit of the boundary layer until the fol- 
lowing section. In the case of the smooth plate, the region of failure of equation 
(1.6) near y = 0 may be identified with the laminar sub-layer, and the failure 
may be explained as due simply to the fact that in this region y is not small in 
comparison with K. For the rough plate matters are more complicated; 
however, we can certainly state that, in this case, the putative origins for the - 
velocity profile, the shearing stress and the mixing length will, in general, all 
differ, and it is very unlikely that equation (1.6) will be admissible in the 
immediate neighborhood of y = 0. 

We may allow for failure of equation (1.6) close to y=0O by relaxing 
condition (2.4) and assuming that 


40, RS gd =e er (2.7) 


where yy is some small value of y, which remains free to vary with x. To preserve 
the similarity which permits use of equations (2.2) and (2.3), it is necessary that 
Yo be of the form yy = y, x/@-2™), with y a constant. 

Such a relationship implies, roughly, that the laminar sub-layer [or, in the 
case of the rough plate, the region in which equation (1.6) fails]-occupies a 
constant fraction of the thickness of the boundary layer. It would be more 
realistic to expect that the absolute rather than the relative thickness of this 
region should remain constant‘). Nevertheless, despite the slight inconsistency 
in this boundary condition, we are able, with its use, to proceed to a similarity 
solution which retains most aspects of the physical problem, and apparently 
gives good results for all but the innermost zone of the boundary layer. 

4) The well-known result that y us/v 10 at the outer limit of the laminar sub-layer implies 


pi Palas thickness varies as w;1. It follows that, for m = 7/8, sub-layer thickness increases with 
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3. Convergence of the Shear Integral 


Integrating equation (2.3) gives 


= 
S 


I 
Oo 
| 
bh) 
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ae 
ee 
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Using this in equation (2.2) and integrating, we obtain 


x 


1 
or ‘| E xX =) 1 “| ay = G he al ' (3.2) 
x 


0 


We introduce t, by means of equation (3.3) 


fi ae (3.3) 
Evidently 


a (=) a (2 -2m)/(3—2m) = of a(2—2m)/(3 = 2m) : (3.4) 


where uw, and c; have their usual significance of ‘friction velocity’ and ‘local 
drag coefficient’. It follows from equations (1.3), (1.4), (1.5), and (3.3) that 


= {pgm Gn? 

ty = (2 7) es : (3.5) 
T, is the analogue of shear stress in our system of reduced variables. Denot- 

ing by Tt,» the value of t, at 7 = mp (i. e. the reduced wall shear), we have 


1 x 
1 ig : . 
isa if E x = H “| dy. (3.6) 
i) 6 


Obviously t, 9 is finite in the physical application, and it is accordingly 
necessary that the integral on the right hand side of equation (3.6) be con- 
vergent. We designate this integral the ‘shear integral’. It seems certain that 
this integral is divergent for values of m of interest in the present context, 7/ 
we retain (2.6) as the outer boundary condition governing the motion. The following 
short investigation of the divergence of the integral is quite non-rigorous, but 
does give some insight into the problem. 


x 
Firstly, it will be noted that, of 7 y and i n dy, x is the dominant term, 
usr 


and especially so in the region close to y = 1, if 7 >oo as y >1. We may 
therefore reasonably assume that the shear integral is convergent or divergent 


1 
according as iL nx dz is convergent or divergent. We introduce y' to denote 
0 
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(1 — x) and we suppose that as 7 > 1 (i. e. 4’ > 0), 7 > co as does (y’)*. Then 
for v' small, we may replace 7 y by 7, and it is seen that if n dy’ can be convergent 


0 
only if « > —1. Also 4 oo as yx’ > 0 [cf. conditions (2.6)] only if «<0. 
Therefore 0 > a > — 1. 
Provided k remains of constant order of magnitude (an assumption con- 
firmed later), dimensional analysis of equation (3.2) gives the result 


1 


aoe (3.7) 


Hyon” (PP Pea = 


It follows that the integral can be convergent only for m > 2, a range of 
values for which substitutions (2.1) are not applicable (and which would, in 
any case, have little bearing on the physical problem). We therefore conclude 
that the shear integral under consideration is divergent. (Note that a similar 
treatment of the corresponding integral for the laminar boundary layer gives 
the result that 7 oo as does the inverse error function of y’. Comparison 
with the error function then confirms the convergence of the shear integral in 
the laminar case.) 

The explanation of this difficulty centers on condition (2.6). This, 1f adopted, 
would imply that the turbulent regime extends throughout the whole infinite 
region y > 0 and would neglect v as small in comparison with K in this region. 
It is evident that such a formulation would lead to such small values of K at 
large y that w would approach U so slowly that the shear integral would be 
divergent. Clearly the integral would be convergent if the fact that the ‘vis- 
cosity’ cannot be less than y were incorporated into the formulation (cf. preced- 
ing remark on convergence of the laminar momentum integral). 

The most straightforward way of avoiding this impasse is to retain equation 
(1.6), but to restrict its region of application by replacing (2.6) with 


eS Se es U (3.8) 
and adding a further condition that the reduced wall shear, Tx, 1S given. 


This may be expressed in the form 


Ou Ou 2_9 2 
| (x DE. => 0 al dy = Tx.0 U2 Pa oon eh dg e (3.9) 


For Ty, in (3.9) to be independent of x, it is necessary that Vs = ny oleae 
where Hx 1s independent of x. The revised conditions replacing (2.6) are then | 
(in terms of the reduced variables) (3.6) and 


= Ty vem (3.10). 


= en ae 
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Either 7; or tT, 9 may be taken as the given parameter specifying the motion, 
These two quantities are uniquely related (for given m and k). 

In the following sections 4 and 6 we make appeal to the available experi- 
mental data in order to deduce-an appropriate numerical value of m and an 
appropriate k-function. In this way we reduce the general formulation developed 
here to a particular form both amenable to numerical analysis and in consonance 


with experiment. Once m is fixed, ty 9 can be discussed more fully than 
hitherto; this is done in section 5. 


4. Numerical Value of m 


The preceding analysis implies the following relationships: 
(a) Boundary layer thickness, ys, is proportional to «/@-2™); 
(b) local drag coefficient, c;, is proportional to x~@-2™/@G-2m), 


(a) Boundary Layer Thickness 


It is generally accepted that, for the fully developed turbulent boundary 
layer over the flat plate (smooth or rough) at zero incidence with zero pressure 
gradient, y; oc x°*®. Data on this point for smooth plates due to ScHULTz- 
GRuUNOW [3] and to TOWNSEND [4] (virtual origin — 21 cm) are shown in Figure 1. 
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é Ls ‘ Fe) 
Thickness of turbulent boundary layers with zero pressure gradient. Fitted lines are for yg oc #°"*. 


+ TownseEnp [4]; lf Corrstn and KistTLer [5];, W SCHULTZ-GRUNOW [3]. 


Both do, in fact, fit this relationship quite well. The data of CoRRSIN and 
KISTLER [5] for a rough corrugated surface also fits the 0-8 power law fairly 
well, when the origin of y, is taken at the top of the corrugations and the virtual 
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origin of x at — 36 in. However, the confirmation of the 0-8 power law, which 
we provisionally accept here, is not conclusive, since much of it relies on inferred 
virtual origins of x. The 0-8 power law for boundary layer thickness corresponds 
to m = 0°875. 


(b) Local Drag Coefficient 


Reference [6] (p. 439) summarizes most of the data available on drag 
coefficients for smooth flat plates at zero incidence with zero pressure gradient. 
In the range 2:5 x10° << R < 10? (R= U x/), c; is found to be proportional 
to x-%20 to a good approximation. However, for the full range of data 
2:5 «10° << R<5x108 a better fit is obtained when the exponent is made 
— 0:17. 

Primary data on drag coefficients for rough plates are not available to the 
author. However, the Schlichting interpolation formula (reference [7], p. 381; 
[6], p. 449) for the range 10? < x/e < 108, where e is the ‘equivalent sand 
roughness’, 


F x \-2:5 
ch = (2:87 +. 1:58 log = (4.1) 
may be recast into the form, 
c, = (function of ) - x, (4.2) | 


where the best fit is to be found with « = 0-193 (see Figure 2). 

It follows that the most suitable values of m would be (i), for smooth plates 
with 2:5x10®°< R< 10’, m=0-875; (ii), for smooth plates with 
2:5 x10®° << R< 5 x108, m = 0-898; (iii), for rough plates, m = 0-880. 

These various types of data are all satisfied reasonably well by the value 
m = 0-875. Accordingly, we adopt m = 7/8 in the developments which follow. | 
Additional data may lead to some change in the value of m, but this is unlikely 
to affect the general nature of the results to be obtained here. 

Difficulties may arise when rough and smooth regimes co-exist along the 
one (uniformly rough) plate. However, since the thickness of the laminar 
sub-layer varies only slowly with x (approximately as x1, cf. previous footnote), 
this will occur in practice only exceptionally. 


5. The Reduced Wall Shear, Tx,0 


With m now fixed, we may relate ty 9 to the experimental data. The rela- 
tionships we develop in this section are in no sense final, but indicate the order 
of magnitude of t, 4 and give an overall picture of its mode of variation on 


smooth and rough plates. The subsequent developments do not depend on this 
section in any way. 


i 
2 
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(a) Smooth Plates 
We require to fit the available data on c, to a relationship of the form 
ey 
Cc; = Ax ; (5.1) 


In the range 2:5 x 10° < R < 10%, the data ({6], p. 439) are very well fitted 
by the value A = 0-074 (v/U)®®°. When the range is extended to 


B90. 610 eR <t-5 10%, 


A requires to be increased to the value 0-077 (v/U)®®°, This value gives results 
deviating from experiment by a maximum of about 10%, an accuracy sufficient 
for our present purpose. Now, it follows from equations (3.4) and (5.1) that 
A=2T, 9, so- that 
py \0-20 
Tx,9 = 0-0385 (7-) (5.2) 
This is the relationship we provisionally adopt to specify tT, 5 over smooth 
plates. 


> 


{= 
a@=- 0-193 
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Figure 2 


ScHLICHTING’s interpolation formula (4.1), and power law approximations to it. Crosses represent 
points computed from equation (4.1). 


(b) Rough Plates 


Here we shall make use of SCHLICHTING’s interpolation formula (4.1). In 
| Figure 2 the function (2-87 + 1-58 log x/e)~° is plotted against «/e for the 


range 10? < x/e < 108. It is evident that we can replace this function to a 


good approximation by one of the form B(x/e)*. The best tit; fora = = "0-198, 


is shown on the figure, but it is more convenient for our present purpose to 


F 


ye os 


ZAMP X/31 
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adopt « = — 0-20, giving B = 0-0264. This line is also shown, and is seen to 
correspond-quite closely to SCHLICHTING’s function. We therefore rewrite (4-1) 
as 


ct = 0-0264 (Zio (5.3) 
Use of equation (3.4) then reduces this to 
T¢. p= 0013 28 (5.4) 


which we adopt as the provisional relationship specifying T, 9 over rough plates. 


(c) Transition from Smooth to Rough 


According to [6] (Figure 21.10, p. 449), the transition from smooth to rough 
over flat plates is marked by a much less definite depression of the drag 
coefficient than is the case for pipes. The point of departure from the smooth 
regime may therefore be found approximately by equating values of T, 9 
obtained from both equations (5.2) and (5.4), giving 


Ue 
Tae a 200 ) (5.5) 


in fair agreement with SCHLICHTING’S ((6], p. 454) approximate criterion for 
10 


1 
=] 

10 \\Y ie 

ee eee 

& z 2 

s&s 3 re 

=] rad > 

iS = ray 
§ Ss lees 
41n72 ss |S s+ 
19 SS §8| §8 

2) Sa = 
ex @x| 2x 
S [es RS§| BS 
) SS 99) am x 2 
‘SS Sse) SUS 
cS 5 QS 6 
LZ Te Le] FE 

-5 -4 = = 
10 10 10 3 10 2 107 
—_ cm ——> 
U 
Figure 3 


Vaniation of T*,9 with »/U and ¢, as given by equations (5.2) and (5.4). The fine lines are ‘con- 

tours of Tx, 9, and are labelled with the appropriate values of Ty o- Lhe diagonal broken line indicates 

approximate lower limit of completely rough regime. Approximate limit of (smooth) turbulent regime 

is shown for three values of x. To right of these three limits regime would be laminar. This figure 
is partly schematic and is intended primarily for illustrative purposes. 
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“admissible roughness’ 


’ 


Ue , 
—~= 100. (5.6) 


The dependence of t, 9 on »/U in the smooth regime and on ¢ in the rough 
regime, as predicted by equations (5.2) and (5.4), is shown in Figure 3. The 
dotted limit of “complete roughness’ corresponds to U é/y = 1500, in rough 
agreement with [6] (Figure 21.10, p. 449). Note that the region corresponding 
to large values of »/U is turbulent only when ¥, (and therefore x) assumes very 
large values. A suitable criterion (for smooth surfaces) is that for x < 105 v/U 
the flow is laminar (reference [6], p. 330). 


6. The R-Function 
Rearranging equation (3.2) 


J Xx -fr iz) axl i 


dx 
3 ees 2 i 1/2 ,.m 
( m)? 9 ae 


hem (6.1) 


It follows that experimental data on the velocity profiles over flat plates 
(zero pressure gradient), which reduce to a relationship between 7 and y for 
each profile, can be used to deduce the k-function. 


x 
The term f y dy proves to be of order of magnitude 10% of 7 y. It may be 
0 


replaced for present purposes by (1 — f) 7 y, where f is taken as constant). 
We adopt the value 6 = 0-88, appropriate to boundary layers characterized by 
To Values about 0-005. Then m = 7/8, 8 = 0-88, lead to the following equa- 
tions, which are modified forms of equations (3.2), (3.6) and (6.1), respectively: 


1 
dy \2 
0-704 a. nnd = (hohe), (6.2) 
Xx 


5) This ‘B-approximation’ introduces the possibility of two distinct types of error. The first is 
that resulting from taking B constant for any fixed value of Tx, 9. Check computations comparing 
velocity profiles (i) calculated with constant f and (ii) using the ‘exact’ equations indicate that 
errors of this type are negligibly small for present purposes. 

The second type of error, namely that arising from adopting the same value of B for all values 
of Tx, 9 raises more complicated issues, but again the errors appear to be unimportant in the present 
context. Strictly, 6 tends to decrease slowly as Tx, 9 increases. 

Note that y and 1/T, profiles (cf. Figures 5 and 7), computed for a given 7, are independent 


x 
of the value of 6 adopted. However, matters such as the relative magnitudes of 7 y and / n dy in 
’ F i) 
any particular boundary layer will be studied most reliably by retaining the exact equations. It 


would have been possible to proceed with the exact equations; but the small improvement in 
accuracy would not seem to justify the great increase in the labor of computation. 
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A 
T4,9 = 0°704 f my dy, , (6.3) 
0 
1 1/2 
(/ nY ix) | 
— . x — 
= 0,839 tm (6.4) 
gis —* 
dn 


Since dy/dyj enters equation (6.4), the notoriously inaccurate process of 
differentiating experimentally determined functions must be carried out in 
evaluating k. This appears to set the limit to the accuracy with which k can be 
found. There is also some imprecision in the values of ys; this tends, in par- 
ticular, to make determinations of k uncertain in outer parts of the boundary 
layer. ee 

Values of k as a function of y/ns, computed from (6.4) for various experi- 
mental data, are shown in Figure 4. We discuss these data, and the k(y/n5) 
function we adopt from them, in the following paragraphs. 


a oKDT 


a4 


Figure 4 


Mixing length function k(y/175), computed for m= 7/8 from experimental data. Points K1, S1 

represent outer limit of influence of laminar sub-layers. Point K D1 represents outer limit of inflecaon 

of smooth surface (see text). Solid curve represents k-function (6.5), adopted in the subsequent 

computations. Broken curve represents k-function (6.6). + KLEBANOFF[8]; W ScHULTz-GRUNOW 
[3]; @ Kiesanorr and Dieu [9]. 


(a) Data of Klebanoff 


KLEBANOFF [8] measured the velocity profile across an artificially thickened 
boundary layer in a wind tunnel specially constructed to give a low turbulence 
level (about 0-03%). In view of the care of the measurements, and the (relative) 
accuracy with which the present author could differentiate the fairly large- 
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scale plot of the velocity profile, these data (with those’of KLEBANOFF and 

DieEHL [3]) must be regarded as the most reliable for our present purpose. 
KLEBANOFF deduces c; = 0:00283 from his measurements of: turbulent 

intensity close to the wall; this corresponds to ty» = 0:00476 cm1/®, The value 


. . 
of T,,9 obtained directly from the velocity profile by evaluation of ui n.4 dy 
0 


li. €., using equation (6.3)] is 000481 cm1!>, Such goodagreement between Ty 
values found by the two independent methods, encourages us to pursue the 
implications of the present formulation. For these data, ys = 0-0616 cm, 

The point labelled K1 on Figure 4 (y/7; = 0.011) appears to mark the 
outer limit of the influence of the laminar. sub-layer on the Klebanoff profile. 
k decreases as the sub-layer is approached, and is about 0-1 in most of the 
sub-layer. 


(b) Data of Schultz-Grunow 


SCHULTZ-GRUNOW [3] gives data for a number of stations along a turbulent 
boundary layer on a flat plate (zero pressure gradient). The velocity profiles 
exhibit similarity to a fair accuracy, and, for simplicity, we here use only the 
data for the station at 530cm. Regrettably, ScHULTz-GRUNOW’s data are 
presented in small diagrams, and the present author does not regard the 
deduced values of &, shown in Figure 4, as having the reliability of the other 
plotted data. ; 

SCHULTZ-GRUNOW measured the drag by means of a floating element 
mounted in the plate. It may be inferred from his data that ci at 530 cm was 
0-0254, which leads to ty 9 = 0:00445 cm1/*. Determination of t, .) from the 
velocity profile via equation (6.4) gives 0-00444 cm1/>. Again the two independ- 
ently determined values of T, ) agree very well. For these data, = 0-0557cm""*. 

The two points on Figure 4 indicated by S1 appear to represent the outer 
limit of the influence of the laminar sub-layer for this profile. 


(c) Data of Klebanoff and Diehl 


To test the k(n/y,) relationship discernible in the Klebanoff and the Schultz- 

_ Grunow data, it was desirable to evaluate k for a surface exhibiting somewhat 

larger values of Ty » and 73. One point to be resolved by such a test was whether 

k should be treated as a ‘universal’ function of /n3 or of 7. The previous 

analytical results follow equally well if & were taken as a function of 7 rather 

than of 7/7. Also, for such data to be reasonably homogeneous with the 

Klebanoff and the Schultz-Grunow data, it was desirable that the turbulence 
level of the free stream be small. 

Data available in the work'of KLEBANOFF AND D1eEHL [9] satisfy both con- 

ditions quite well. Their work was carried out in the same tunnel as that of 
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KLEBANOFF, and the data to be used here gave a value of 4; about twice those 
of the previous data. 

KLEBANOFF and DIEHt covered the leading two feet of the test plate with 
No. 16 sand paper, so as to artificially thicken the boundary layer, which would 
be studied in detail sufficiently downstream for the layer to have resumed 
equilibrium characteristics. Our interest here is in the velocity profiles at a 
station 0-08 ft downstream of the roughened surface. The boundary layer over 
the roughened surface was apparently turbulent right from the leading edge, 
and this 2:08 ft profile can be considered characteristic of the roughened surface 
(and in equilibrium with profiles over the roughened surface), except in a sur- 
face layer of thickness of order 0-1 in. 

KLEBANOFF and DIEHL give data for free stream velocities of 35, 55 and 
108 ft/s. All profiles agree closely, and here we have used a mean curve through 
the three sets of data. On the present interpretation of these data, 7; = 0-105 
and T, 9 = 0-0134. The point labelled KD1 appears to lie at the inner limit of 
validity of the assumption that the profile is in equilibrium with the rough 
surface. Note that points within the region of influence of the smooth surface 
would give abnormally high values of k, if they were treated as homogeneous 
with the outer part of the profile. 

It is seen that the data of KLEBANoFF and DIEHL lead to a k(y/n5) function 
which does not differ significantly from that suggested by the data of KLEBA- 
NOFF and SCHULTZ-GRUNOW. It also appears that our assumption that k varies 
with (7/3), rather than with 7, is justified. 


(d) The Adopted k-Function 


Clearly, k is about 0-36 in a region close to the wall, and decreases mono- 
tonically to about 0-14 in the outer parts of the boundary layer. The data are 
scarcely adequate for us to go beyond this. However, to explore numerical 
solutions of the boundary layer equations, we need to specify & in precise form. 
On the basis of the experimental data, the following k-function (for low free- 
stream turbulence) appears appropriate, and is used subsequently in this work. 


6s Zoe pen0ss, | 


Ns 
02<1<06, &=0-0871 (2) 
"6 "6 ; 


(6.5) 
0-6 < sae a1, ke 01355. 


(6.5), which is shown as the full curve of Figure 4, is in a form chosen 
because it was convenient for the numerical work we require to perform in 
integrating equation (6.2). Certain inelegancies in (6.5) can be avoided in other 
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formulations, such as 


OS sre Da, 086" 
6 


(6.6) 


OF tl <r  Fiep es (4) 
Us 5° 6 i 


which fits the data nearly as well. The broken curve of Figure 4 represents (6.6). 

These forms are, of course, tentative. Note that the three and four signif- 
icant figures of (6.5) are not to be taken seriously. Calculations were begun 
with a B-value which was later amended. The numerical values of (6.5) were 
needed so that these calculations could be retained. 


7. Numerical Solutions of Equation (6.2) 


The method of solution of equation (6.2) subject to conditions (3.10) and 
(6.3) is, essentially, to regard 7, as a datum, and to proceed by stepwise inte- 
gration from 7; towards 7. Values of t, 9 and 7, emerge as part of the solution 
obtained. The method is similar to that used in integration of the concentration- 
dependent diffusion equation [10], in that the y-intervals used in computing 


1 
successive values of ff » 4 dz overlap those used in computing the decrements 


x 
in». Near y = 1, 7 = nz, the integration is started analytically by means of the 


approximation 
1 


fa x 44 = (7H H5)"? (1 — x), (7.1) 


x 
which leads to equation (7.2) for k = 0-1355 


Pl = yl — 0-1211 x5 /* (4 — x)”. (7.2) 


Solutions of equation (6.2) were obtained for the following values of 7s: 
0-05, 0-06, 0-08, 0-10, 0-12. These give a range of values of Ty \9 which includes 
the available experimental data, and which can be expected to contain most 
cases of practical interest. The number of steps dividing the x interval, 0 to i, 
varied from 41 to 58. 

The primary numerical results obtained consisted of the 7(z) functions, 
_ which are shown in Figure 5, and the set of values of 


1 
[unde 
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Figure 5 


Experimental and computed velocity profiles in reduced form 4(y). Figures on computed curves 
denote appropriate values of Ty, 9. W ScHULTZ-GRUNOW [8]; + KLEBANOFF [8]; M CorRsIN and 
KisTLER [5]; @ KLEBANOoFF and DIEHL [9]. 
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Figure 6 


The relationship between Ty, and 76. Crosses denote computed points. Curve represents equation 
(8.1). 
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Figure 7 
Computed shear stress profiles for the five values of Tx, 9 Shown. 


Figure 8 
Computed shear stress profile for Tx, 9 = 0:00489, compared with experimental points for 
Tx, 9 = 0:00476. + KLEBANOFF [8]. 
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which occur in the course of the computation. Also, from equations (3.5) and 
(6.2), 


i 
1, = 0-704 ih nx aX (7.3) 
x 


so that the profile of t, follows at once. Figure 6 shows the Ty ,9 (ys) relationship 
found in this way, and the various T/T» (or T/T) profiles are given in Figure i 


8. The Various: Properties of the Boundary Layer 


Other boundary layer properties may be computed from these primary 
numerical results: The various properties of the layer, as computed, are pre- 
sented and discussed briefly here. 


(a) Shear Stress Profiles 


The shear stress profiles':(Figure 7):.are all-sigmoidal, but not identical. 
KLEBANOFF [8] measured the: turbulent: shearing stress in a boundary layer for 
which Ty 9 was:0-00476. Figure 8 compares his data with the computed profile 
for Ty 9 = 0:00489; the ys adopted for KreBanorr’s data was 3-12 in. The 
agreement between analysis ‘and experiment is seen'to be fairly good, and is 
rather better than that obtained by CuausER’s analysis ([{11], Figure 16). 

In Figure 6 the computed points are: compared. with the parabola passing 
through the point computed for; = 0-08, and representing the relationship 


Tyo = 1:352 7p. (8.1) 


(b) Velocity Profiles 


Figure 5 shows various experimental data, as well as the computed reduced 
velocity profiles. The profiles of KLEBANOoFF and ScHULTz-GRUNOW, used 
previously in establishing k, are given. The 55 ft/s profile of KLEBANOFF and 
DIEHL is shown, the other similar ones being omitted for clarity. Also shown is 
the profile of Corrstn and KIsTLER [5] at station 102” of their boundary layer 
along a rough corrugated surface. The agreement between observation and 
analysis is especially notable in the case of the Corrsin and Kistler data, since 
these did not enter into the establishment of &. 

The computed profiles agree about as well as do observed profiles with the 
velocity profile ‘laws’ proposed from time to time. 

(i) Power law velocity profile. Plots of logy) y against logy 4/ns for each 
computed velocity profile show that an appreciable region of the layer might 
be supposed to obey a power law. The exponent increases practically linearly 
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with tT, 9 from 1/9-8 for the ‘smoothest’ profile to 1/1-7 for the ‘roughest’. The 
exponent corresponding to much of the smooth plate work (T%.9 © 0:0046) is 
about 1/7. ; 

(ii) Velocity defect law. Plots of (U — u)/u, against a reduced abscissa 
related to distance from the wall, such as y/y,, have been found [11-15] to give 
groupings of points which might suggest the existence of a ‘universal’ defect 
law holding in an outer region of turbulent boundary layers along both smooth 
and rough plates. 

Such formulations seem to imply two different types of universality. These 
are, firstly, a universality embracing profiles occurring at different stations 
along the one boundary layer, and, secondly, an universality embracing profiles 
occurring on different boundary layers, i. e. on smooth plates characterized by 
different values of (y/U) and on plates of various roughnesses. 

Universality of the first type cannot be exactly consistent with the present 
formulation (m + 1), which implies similarity in velocity profiles along the one 
layer. (The data of ScHULTZ-GRUNOW [3] is, perhaps, the best available for 
comparing velocity profiles at a number of stations along a plate at zero inci- 
dence. Curiously enough, these data are often supposed to illustrate the velocity 
defect law. Yet it can be shown, at least as persuasively, that his profiles exhibit 
similarity to a good approximation.) Universality of the second type, however, 
can be examined further by means of the present approach. The corresponding 
plot here is of (1 — )/t1/*, against »/7,. In view of equation (8.1), we use here 
the equivalent, and more convenient, (1 — y)/n, against 7/n;. Plots of the 
computed profiles are given to a linear scale in Figure 9, and with n/n; loga- 
rithmic in Figure 10. 

Figure 9 should be compared with, Figure 10.5 of [15] (p. 245). The spread 
of the experimental points in TOWNSEND’s figure is about as great as the 
deviations between computed curves for the appropriate t, ) values. Admit- 
tedly the x-values for the experimental data are not all equal, but they tend 
to be of the same order of magnitude. There are also factors tending to obscure 
the spread of points at equal x; for example, the profiles for CorRsin and 
KIsTLER’s rough surface were at x-values very much less than those appropriate 
to the smooth surface data of SCHUBAUER and KLEBANOFF [16]. Similar remarks 
hold for Figure 10, which may be compared to Figure 10.4 of (15] (p. 244) and 
Figure 5 of [11]. 

(iii) Logarithmic profiles and roughness effects. Another interesting aspect 
of Figure 10 is that, at any fixed t, 5, the velocity is a linear function of logy/ns 
to a good approximation in the region 7 < 0:2. The present analysis thus 
connects with the proposal [11, 13, 15] that this part of the velocity profile is 
logarithmic. 

Also, the various curves are roughly parallel in this ‘logarithmic’ region. 
This is at least consonant with the concept [12, 17] that roughness changes 
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Figure 9 


Velocity defect law, usual formulation. Linear plot of computed velocity profiles for the five values 
of Tx, 0 shown. 
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Figure 10_ 


Velocity defect law, usual formulation. Logarithmic plot of computed velocity profiles for the five 
values of Tx, 9 Shown. 
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the velocity profile in the ‘logarithmic’ region by decreasing u/u, by a constant 
amount. : 

(iv) An alternative velocity defect law. It is simply shown that use of 
approximation (7.1) in equation (6.2) leads to the result that velocity profiles 
may be put into the form 

fad: 
ent 


bac . 1) 
= function (z-) ; (8.2) 

Although agreement between plots of the various profiles depends solely on 
the efficiency of approximation (7.1), which is accurate only in outer parts of 
the layer, this form of defect law is found to work quite well through the whole 
depth of the layer. 

A particular form of (8.2) appropriate to the outer part of the layer may be 
found by rearranging equation (7.2), using the binomial theorem, and neg- 
lecting terms in (1 — »/n 5) of higher order than the first. The result is equation 
(8.3): . 


L fa nH \2 
aie Sectgusg ( = 7 oa (8.3) 
This is similar in form to HAma’s [17] empirical square law for the outer part 


of the boundary layer. It supports HAma’s second power rather than Ross’s 
[18] 3/2 power. 


(c) Profiles of Eddy Viscosity 
Equations (1.4), (1.5) and (2.1) yield 


ee 2 oe ay rr (2m—1)/(3—2m) 
K = k Y] an U x (8.4) 
: so that we may introduce a reduced eddy viscosity K, by the relationships 
K aC BY (2m —~1)/(3 —2 m) ) pao! 2.2m dy : (8 5) 
= KK, UO x , Aves he a ae : 


It follows from equation (3.5) that 


Higeenen Ps (8.6) 


which enables profiles of K, to be computed simply. These are shown in 
4 Figure 11. They are fairly similar in shape, the scale of K, being roughly 
_ proportional to ty 9. Since k? enters the second of equations (8.5), one expects 
the discontinuities in the adopted k-function at /73 = 0-2 and 0-6 to be reflected 
in the K, curves. The slight dip in the curves near n|ns = 0-6 is explained in 
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this way, and it is likely that the-peaks near ning = 0:2.are rather too sharp for 
the same reason. 


0 0:02 0:04 0-06 0:08 0:10 OT2 
7] cm? —— 


Figure 11 
Profiles of eddy viscosity in reduced form K,(y), computed for the five values of Tx, 9 shown. 


(d) Velocity Component Normal to the Plate 


Profiles of reduced velocity normal to the plate are simply computed, using 
equation (3-1). The profiles are found to be typically sigmoidal, but, as one 
would expect, are not exactly similar to each other. 


(e) Relation Between Roughness and Drag 


The numerical solutions also yield a relationship between 7 and Ty 9. At 
the ‘rough’ end of the t, 9 range, 7) varies approximately as t% 9. Taken 
with equation (5.4), this suggests the reasonable result that, for the rough 
regime, 7) may be proportional to «. 


9. Free Stream Turbulence and Pipe Flow 


The k-function adopted here is intended to apply to the case where the 
turbulence level in the free stream is small. As free stream turbulence increases 
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one would expect & to increase, particularly inthe outer regions of the boundary 
layer. Turning to the question of turbulent flow in pipes, it may be reasonable 
to regard the pipe radius as the boundary layer thickness, and the equivalent 
‘free stream turbulence’ as some representative level of turbulence along the 
opposite radius of the pipe. Thus the k-function for pipes might be expected 
to be rather like that for a boundary layer with a high level of free stream 
turbulence. This is, in fact, what is found for the mixing length data of NIKU- 
RADSE [19], when this is recast into the form k(y/y,) = 1-78, 


10. Discussion 
(a) Universal Functions and Occam’s Razor 


WILLIAM OF OccAM pronounced, ‘entia non sunt multiplicanda praeter 
necessitatem’. Fluid mechanicians may translate: ‘Universal functions must 
not be multiplied beyond necessity.’ For each, separate universal function, as 
used in fluid mechanics, is simply a separate hypothesis about the phenomenon. 

In the present work, we have reduced description of the turbulent boundary 
layer, in all its gross particulars, to a single ‘mixing-length’ function®). The 
alternative formulations of CLAUSER, CoLEs and TOWNSEND apparently 
require one universal function to describe the outer velocity profile, and 
another to describe the inner velocity profile in the case of smooth surfaces 
(law of the wall). CLAUSER and Hama have made some suggestions as to how 
to cope with the problem of the inner part of the velocity profile over rough 
surfaces, and this, in effect, involves yet a third hypothesis. The present 
formulation, with its single hypothesis, is relevant to both smooth and rough 
surfaces, and involves no problem of joining inner and outer profiles. Note, 
also, that no further universal function is required to generalize the present 
approach to the case of non-zero pressure gradient. 


(b) Downstream Variation of Drag and Layer Thickness 


A further aspect in which alternative formulations appear to be inadequate 
is in the manner in which they fail to take full cognizance of the downstream 
variation of local drag coefficient and boundary layer thickness. Thus in [11] 
é (p. 18) a universal defect law is made the starting point of an analysis of the 
boundary layer. Immediately a difficulty arises because u, varies with x, but, 
because the variation of u, with x is relatively slow, the theory is developed on 
the basis of u, constant, i.e. ‘pipe-like’ elements are introduced into the 
analysis at the outset. 


8) We do not, of course, pretend to give an adequate description of the laminar sub-layer, 
where it exists. However, this seems a small price to pay for the advantages of the present for- 
_ mulation. 
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Another form of this tendency to treat boundary layers as pipes is inade- 
quate recognition of the loss of momentum suffered by each parcel of fluid 
during its passage along a thickening boundary layer. The terms on the left 
hand side of equation (1.2) relate to the rate of this loss of momentum, and 
clearly form a most important element in this equation. They do not,.of course, 
arise in pipe flow. 

It is quite obvious that these terms are greatly affected by the rate of 
growth of boundary layer thickness, and it follows that the profiles of velocity 
and shear stress will be quite different if, say, y; x x°* instead of y; ox x°°. 

CLAUSER [11] attempts to represent the outer regions of constant pressure 
layers by velocity profiles of laminar boundary layers with various non-zero 
velocities at the wall. He wishes to reproduce the effect of constant eddy 
viscosity in the regions, but his results are appropriate to a layer growing in 
thickness as x°5, and having momentum loss characteristics of a laminar and not 
a turbulent layer. His conclusion that the viscosity in the outer 80% to 90% 
of the layer is constant and proportional to the product of U and the displace- 
ment thickness (i.e. to U x°8, with m= 7/8) is, therefore, not valid. The 
present analysis predicts K «x U x8 and gives profiles of K rather more in 
keeping with experimental data [4}. 


(c) Physical Implications 


The physical implications of this work will be explored further in a proposed 
second paper giving the extension to non-zero pressure gradient. Here it suf- 
fices to state that the momentum considerations introduced in equation (1.2)§ 
together with the adopted mixing-length hypothesis, are capable of predicting 
and connecting the gross characteristics of the turbulent boundary layer. 

Objections to the mixing length approach, such as those of TowNSEND [4], 
have substance, of course, and are a useful corrective to a literal belief in the 
physical existence of the simple mixing-length model. Yet, although the model 
cannot be exact, it 7s an approximation to mixing processes which do take 
place, and it is not unreasonable to expect that a mixing-length function might 
be found which describes gross features of the motion to a fair enough accuracy. 
This is apparently what we have done here. | ; 

The full scope of the present approach has not yet been explored. Apart 
from extension to non-zero pressure gradients (which shows, inter alia, that 
the k-function used here leads to results in agreement with observation of 
turbulent layers in positive and negative pressure gradients, including layers 
at separation), there remains the question of flows in pipes and closed and open 


conduits. There is some evidence that the approach may prove valid in this 
somewhat wider context. . | 
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Zusammenfassung 


Das hier entwickelte, formale, mathematische System beschreibt turbulente 
Grenzschichten an glatten und rauhen Platten und vermittelt eine zusammen- 
hangende und folgerichtige Behandlung von solch groben Schichteigenschaften wie 
Durchschnittsgeschwindigkeitsprofil, Schubspannung, Wirbelviskositat und Grenz- 
schichtdicke. Die erhaltenen Lésungen zeigen gute Ubereinstimmung mit den Ver- 
suchsergebnissen. Die gegebene Formulierung kann in gewisser Hinsicht als eine 
neuartige Mischungsweghypothese betrachtet werden. 

Die Arbeit beschrankt sich auf Grenzschichten mit Nulldruckgradienten. Eine 
Verallgemeinerung auf von Null abweichende Druckgradienten ist einer spateren 
Mitteilung vorbehalten, in welcher gezeigt werden soll, dass die benutzte Methode 
zu mit dem Experiment iibereinstimmenden Ergebnisse fiihrt und keiner weiteren 
Hilfshypothesen bedarf. 

Die physikalische Bedeutung der Arbeit kann darin gesehen werden, dass ein 


1 allerdings nicht ganz genaues Mischungswegmodell dennoch auftretende Mischungs- 


vorginge naherungsweise beschreiben kann, was somit die Aufstellung einer 


Mischungsweghypothese erméglicht, die mit geniigender Genauigkeit ein grobes 
- Bild verschiedener Aspekte der beobachteten Vorgange gibt. Es ist somit méglich, 


zu einer handlichen und folgerichtigen mathematischen Formulierung zu gelangen, 


- mit deren Hilfe allgemeine Eigenschaften von turbulenten Grenzschichten an festen 


Oberflachen untersucht werden konnen. 


_ (Received: Dezember 8, 1958.) 
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Resistance on a Circular Cylinder Due to any Number 
of Vortices Lying in two Rows 
By Miss Durca Roy, Calcutta, India*) 


‘ Introduction 


Consider a circular cylinder in a uniform flow with a wake of vortices. The 
purpose of this paper is to calculate the resistance experienced by the cylinder. 

First, a general formula is obtained with no assumptions about the vortices, 
except that half of them are of strength k, and half of strength —k; the result 
is contained in equations (11) and (16), under the further assumption that the — 
vortices lie symmetrically in two lines parallel to the uniform flow. The formula 
for the resistance is made more explicit in the double series shown in (20) 
and (21). 

The present work differs from what has been done previously in the fact — 
that it yields an exact expression for the instantaneous resistance, instead of 
merely a time-average. This is, however, achieved only by limiting attention 
to the case of a circular cylinder. With reference to calculations of the time- 
average, the usual treatment [2]*) involves the replacement of a semi-infinite 
wake by a fictitious wake extending to infinity in both directions, but this 
procedure may be avoided {1}. 


Description of the Present Work 


Where there are two rows of vortex filaments, in one of which each filament 
is of vorticity + and in the other of vorticity —k, placed symmetrically — 
to one another; and where a cylinder of radius C has been introduced with its — 
centre at the point of origin in a uniform flow along the x-axis, the complex 
potential is found with the help of the circle theorem and is given by 
tk 
an 


x [y log (z — z,,) — D/ log (z — Zy) — Dy log (— — 2m) . Ds log(< — zy)] ; 


: et bes at the School of Theoretical Physics, Dublin Institute for Advanced Studies, Dublin, 
reland. : 


2) Numbers in brackets refer to References, page 508. 


W= -U (2+ =) + 


eS ae 
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where 


2m = Km FTV ANd By = ty — iy. 


W can be written as 


W—=—U (z+ 5) 454% 


/ (1) 
>: [log (z — z,,) — log(z — z,) — log(z - =) + log(z — :-) + log 2 . 
m / ‘pr mn 


If X, Y denote the components of thrust acting on the cylinder along the axes 
x and y then, 


X-iY=-ifp paz, (2) 


where integration is along the boundary of the circle z z= c?, where c is the 
radius of the cylinder. 
The pressure # is obtained from BERNOULLI’s equation 


fi) 1 
f= -Set Fl), (3) 


where gq? is the square of the velocity, y the velocity potential and g = —V 9. 
To evaluate F(é) we see that at infinity = const = JT say, 


g=WU ‘and SP _ 0 for all t. 


Therefore (3) becomes 


1 1 ay : 
p=H-=eU-— Feet (4) 
Now, 
dw \ [ dw Op ow Ow 
ge = (%) (%), A some, te ee 6) 


Substituting the value of p from (4) with the help of (5), equation (2) becomes 
: Higgs Ow \2 ‘ Op\ ( Oz 


giving the components of the total force on the cylinder, The result (6) is 
obviously an extension of the Blasius result for the steady case, 
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Ie Ow \2 
D Ges Ph Gs Sa | (s-] Oz (6a) 
and ‘ ie 
F . a 
Kosa A | (SP) (52). (6b) 
Therefore . 
| EE Ba) ee A Oe. (6c) 


It is easy to identify X,— 1 Y, as the complex thrust on the cylinder due to 
two rows of vortices when the vortices are stationary. 

On the other hand, X,— 7 Y, is the contribution to the complex-thrust 
from the unsteady character of the flow created by the movement of the 
vortices. 

From (1) we get 


0 tk 1 HE 1 1 
=a pa Piers mem rst aie 


where 
= Ce O : 
Canoe a and a, = ee (7) 
With the help of (7) equation (6a) becomes 
igi 1 1 1 1 
4-3 Y= 2m PIP eee ee | — hy hare (8) 


viexve(S-32 7 =): 


where m, denote the vortices on the upper row and #’, g’ on the lower row. 
Simplifying (8) finally we get 


zm a Zim Zp.) aCe (an =i Bpi) 
xy (Tz oe Fm Zp. — ©) 


‘BE a +o) — ty te By + Fa) ©) 


=F = Za & (ig c?) 


nis ES ea 
eipa ore 


™m ™ p ! 


ee ee ee 
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If the vortices on the upper and lower row move symmetrically, then (8) can 
be written as 


3 ; k? 
X,-i¥Y,=2— yy E a —- aie ee ca as ig 


=m — %&p Em — %&n, & 


{ 1 1 
Vem — tn) Em = On) Em — On) Em — By) i 
+iokUC(Y = - 3S a) 


As the vortices are symmetric in (10), replacing m by v in the first exponent 
of (10), adding and dividing by 2, and simplifying, we get 


(11) 
{ 1 1 \ 
3 VG en — c*) (Z Zn G7) (Zm an 63) (Z en c*) 
+igKUCY({-=). 


From (1) again, with the help of the boundary condition we get 
dw ik — Em mek en) 
ee Ply e+e +Y zie _| (hy = (12) 


From (12) 


op 1 / ow ee 3 
ot 2 ( Be Tae >) 
and from (6b) and (13) we get 
Z : ea OO Oz 
X,-i¥,=—ie/ (+7) Ce) a, (14) 
Eliminating z from the above integral (14) with the help of z z = c® we get 
Zo 
X,—-1Y,= -iKeor(y fe — ¥ 2). (15) 


Adding (9) and (15) we finally get X —7 Y. 
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In the case when the vortices move symmetrically then (15) becomes 


Ds Za 
X,-i%=ikeot (Y Zz Saal (16) 


Summing (16) and (11) we finally get X —7 Y. 
In order to find out the velocities of the vortices we take the classical 
velocity which is given by 


Ow, 


Oz eg a lm — 8 Uin =n — 9.9 = Lay (17) 


where 
tk 
i a, log = 2) . 


From (1), when the vortices move symmetrically, with the help of (17) we get 


= ea: 1 1 z, Zz, 
LA ieee i) Zm —% Zim —Z, es aoe ee. a 


m an 


+ Ue 


and 
: Wh eee i a z, Z, 
=m ab Ee a? er z, er a oe Sata] 
(18) 
c 
+uU(1- =) 


EL OR Sh CP Aeiy Spi Sey Ay eee eee. (2m + Fm) (19) & 


mn 
m nN Zin Zin (c? — 2m, Zn) (c? Aeon Zn) 


— YS Gin Stn Fin) = Fm Em En En (Em + Fm) 
m Nn 2m =m eo eee) (C? — Bn Zy) 
a 1 1 
iKeguC Y (= =) 
m m 


Then finally X — 7 Y is obtained by summing (11) and (19). 


—— ee ee ee 
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Now if the vortices be arranged in such a way that 
2 = Mat > a rk and _ana—+ ib, 


then substituting this value expression (11) becomes 


ae 7 = — nn) (mn a® + B) wn. 


mn a+ B)*+ a? bt (nm — m)*t {(mm a2 + 8)? + as + orc] | (20) (20) 
| 


NE a ee eae 
m* a == 


The expression (2) is convergent, as 


Sy. Z ; (m + n) (mn a* + B) 
m nm er ee 


We ay Ve m+n 
Saher ee Hier 23 me n3 a& ad m3 n> ’ 
, L 1 
a = (are 5 0 erat? 
m n 
But both the series 
i 1 

eds m? n8 = yah pre ge 

m n m n 
are convergent; again, as the series 


if 
ease 


m 


is convergent the U depending series in X,— 1 Y, is also convergent and 
therefore X, — i Y, is convergent, when m and both take an infinite number 
of values. 

Taking the values z,,, z,, etc., expression (19) becomes 


1 1h 
men (m2 at + 7 08) 2ma— (mea? ——b') 2na—mad® 


apes 
m n (m — n)? a? (m= at + 7b] 
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(me at + 408) 2 ma — (m2 at — 7-68) 2 a + mad? 


fo ol oe eae - 
= (om? att = 0°) {(m — n)? a2 + b} 


—_— 


1 
ct{ana(me at — 50) + mad'| — (me ar + 7) (n? a? +, 08) 2ma 


eae cr : : 


2 pe apet eal ie Lr a ak iss 
(m? at + | 0°) e+ (m? a + | 04) (nt at + 704) She (mn a? + 0°) 


ot |2na (ma 7b) — mari — (m2 a + 7) (na? + 7 04) 2ma 


1 =| Lp 2 1 ai ; *( ret ES 3 
(m? a? + 0+) {ot + (meat + 7 8) (n° +70) 2c2|mna Z b 


m 


(mm? a+ 7B) 
4 


~29KUCabd y 


It can also be shown that expression (21) is convergent, when both m and n 
take an infinite number of values, and expressions (20) and (21) are both real. 

Therefore, when there are an infinite number of vortices which move sym- 
metrically on two parallel rows, where the distance between the two vortices 
of same strength on one row is 2 a, and the distance between the two parallel — 
rows is b, then the force experienced by the cylinder in a uniform flow is along © 
the x-axis and is given by X — 17 Y, which is nothing but the sum of the two 
expressions (20) and (21) respectively. 

In another paper the author will discuss the result that is here obtained and 


also the case when the vortices are placed asymmetrically following the laws 
of the Karman vortex street. 
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. . . ; 
In dieser Arbeit wird der unterschiedliche Widerstand an einem Kreiszylinder — 
berechnet, der in einem steten Strom angebracht ist und von . 
Wirbelbewegungen in zwei Reihen begleitet ist. 
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Uber die Kondensation verschiedener Dampfe bei 
adiabatischer Expansion 


Von SERGIO BERNascont') und Grore Gorscu?), Locarno-Monti 


1. Einleitung und Problemstellung 


Verwendete Symbole 
dp Druckanderung; 
I ander Photozelle wirksame Strahlungsintensitat zu einem beliebigen Zeitpunkt 
[erg/cm? - s]; 
I, wie I, jedoch zu Beginn der Expansion [erg/cm? : s]; 
a van’t Hoffscher Faktor; 
M,; Molekulargewicht des Stoffes 7 [g/Mol]; 
m; Menge des Stoffes 7 [g]; 
Py) Atmospharendruck [dyn/cm?]; 
P* bo + 4p [dyn/cm*); 
y Radius der Trépfchen [cm]; 
7,» Kkritischer Gleichgewichtsradius [cm]; 
Ubersattigung; 
kritische Gleichgewichtsiibersattigung ; 
Temperatur [°K]; 
eo; Dichte des Stoffes 7 [g/cm]; 
o,/; Oberflachenspannung von 7 gegen 7 [dyn/cm]. 


In dieser Arbeit wird die Kondensation verschiedener Dampfe an Kondensa- 
tionskernen bei adiabatischer Expansion untersucht. Die in der Literatur bisweilen 
vertretene Auffassung, dass das Wachstum der Kondensationskerne entlang den 
Gleichgewichtskurven fiir Lésungs- bzw. Mischkerne verlauft, geht auf Vorstel- 
lungen von THomson [1]’), KOHLER[2] und JuNGE [3] zuriick. Thre analytische 
Form wird durch die Gleichungen (1) und (2) wiedergegeben, und zwar: 

Fiir Kerne léslicher Substanz durch die Kohler-Gleichung 


oa ee (1) 


mit 
Zi Sv/a My.0 


21,0 R* T *: 


1) Osservatorio Ticinese della Centrale Meteorologica Svizzera. 
2) Forschungsstelle am Osservatorio Ticinese der Eidgendssischen Kommission zum Studium 


der Hagelbildung und der Hagelabwehr. 29 ; 
3) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 518. 
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3M,;,5™m 
Ca H,0 * fiir ideale Lésungen , 
4% Qo Ms 
3 Myo Ms hae 2 
Co—4 2 fiir dissoziierende Losungen . 
42 04,0 M s 


Fiir Mischkerne, das heisst Kerne, welche zum Teil aus léslicher, zum Teil aus 
nichtléslicher Substanz bestehen, durch die Junge-Gleichung: 


Crm is Cy 


ae 
Y vs — xi 


as (2) 


7, = Radius des spharisch angenommenen unléslichen Bestandteils. 
Wie man leicht erkennt, geht Gleichung (2) mit C, = 0 in die Thomson-Gleichung 


Cy 
S= Pa (3) 
fiir reine Wassertropfchen iiber, ferner mit 73 < rv? in die Kohler-Gleichung. Unter 
(771 Sy) Wollen wir den Punkt auf den durch die Gleichungen (1) und (2) gegebenen 
Kurven verstehen, fiir welchen die Bedingung 


dS aS 
Pw eras (4) 


r f$—> 


100% 


Figur 1 
Verlauf der Kohlerschen Gleichgewichtskurven fiir zwei Kerne schematisch gezeichnet. 


erfiillt ist und mit WigELaNnp [4] den Teil der Kurve, fiir welchen 7 < 1,,. ist, das 
Stadium der Vorkondensation nennen, jenen Kurventeil mit 7 > 7 das Stadium 
der eigentlichen Kondensation, Entscheidend dafiir, ob ein Kern Ae das Stadium — 
der ee eintritt, ist nach dieser Betrachtungsweise die Bedingung, dass — 
die bersattigung Ss in der Umgebung eines Kerns mindestens den Betrag der 
kritischen Ubersattigung Sy, erreicht haben muss, wobei diese nur von der Grésse 


——- 


she “step weer ~~ 
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und der stofflichen Zusammensetzung des Kerns abhangt. Die skizzierte Theorie 
beriicksichtigt weder die Diffusion noch die Warmeleitung. 

Nach der gegebenen Darstellung kénnen bei einem Vorgang, durch welchen 
eine kernhaltige Luftprobe um einen-gewahlten Betrag S, tibersattigt wird, nur 
jene Kerne in das Stadium der Kondensation eintreten, deren kritische Ubersatti- 
gung S,, SS, ist. JUNGE [5] folgert aus eigenen Beobachtungen im wesentlichen 
dasselbe Verhalten der Kondensationskerne. Er iibersattigt eine kernhaltige Luft- 
probe durch eine kleine, adiabatische Expansion, fotografiert die dabei entstehenden 
Trépichen und komprimiert sofort wieder auf das Ausgangsvolumen. Diesen Vor- 
gang wiederholt er mit stets grésser werdenden Expansionsstufen. Um ein Aus- 
regnen der Tropfchen zu verhindern, lasst er die Kompression in kiirzester Zeit auf 
die Expansion folgen. Aus den Ergebnissen dieses Experimentes folgert JUNGE, 
dass eine 10prozentige Ubersattigung erforderlich ist, um etwa 50% aller « Kerne 
zu aktivieren» (diesen Ausdruck verwendet man jiiblicherweise fiir das Uberschreiten 
des kritischen Radius). Hingegen ergab ein von WIELAND [4] mit einer Mischwol- 
kenkammer durchgefiihrtes Experiment, dass bei einem natiirlichen Aerosol bereits 
Ubersattigungen von <1,2% geniigen, um 1000 bis 10000 Trépfchen/cm® aus- 
regnen zu lassen. 

Kerndichten gleichey Gréssenordnung wurden am Messort auch anhand von 
langeren Messreihen mit dem kleinen Scholzschen Kernzahler und dem automati- 
schen Kernzahler nach VERZAR gefunden. Die mit einem Scholzschen Zahler und 
gleichzeitig mit einer Mischwolkenkammer ausgefiihrten Messungen von GortscH [6]. 
liessen ebenfalls feststellen, dass weniger als 1% Ubersattigung ausreicht, um mehr 
als 50% aller Kerne zu aktivieren. Diese Ubersattigung ist demnach um eine 
Zehnerpotenz kleiner als jene von JUNGE, welcher Widerspruch einen der Verfasser 
[6] veranlasste, erneut auf die Bedingungen hinzuweisen, unter denen die Gleichun- 
gen (1) bis (4) zulassig sind. Wesentlich an ihnen ist, dass sie Gleichgewichtsbezie- 
hungen darstellen, die die Zeit nicht enthalten, welcher Umstand notwendigerweise 
bei der Deutung von Experimenten beriicksichtigt werden muss. Unterhalb der 
Sattigungsfeuchte kann dieser Vorschrift ohne weiteres nachgekommen werden. 
Tatsachlich sind die Gleichungen (1) bis (4) bis zu Feuchten von nahezu 100% 
experimentell bestatigt worden (siehe ausfiihrliche Literatur unter [7)). 

Bei Uberschreiten der Feuchtesattigung ist die den Kernen zur Kondensation 
zur Verfiigung stehende Zeit beschrankt. Dauer und Ausmass der Kondensation 
werden dabei durch Diffusionsvorgange bestimmt, nach deren Ablauf sich wieder . 
ein Gleichgewichtszustand einstellt. In neuerer Zeit ist verschiedentlich versucht 
worden, dieser Tatsache gerecht zu werden, indem man bei Berechnung der Kon- 
densationsvorgange die Diffusion und die Warmeleitung beriicksichtigte [8, 9]. 

Die diesen Rechnungen zugrunde gelegten Modelle eignen sich jedoch nicht zum 
Vergleich mit den nachfolgend beschriebenen Experimenten, weil die dabei ange- 
nommenen Zustandsdinderungen so langsam verlaufen, dass man sie mit einer 
kleinen Expansionskammer nicht nachahmen kann. " 

Bei dem eben erwahnten Experiment von JUNGE ist die Ubersdttigung bei allen 
héheren Expansionsstufen stets iiber die Aktivierungssattigung der wirksamsten 
_Kerne hinaus gesteigert worden (sofern verschiedene aktive Kerne existieren), und 
gwar in so kurzer Zeit, dass die Trépfchen noch nicht aussedimentieren konnten. Wie 

man aber bei allen Experimenten mit Expansionskammern bis anhin einwandfrei 

beobachtet hat, beteiligen sich unabhangig von ihrer Aktivitat alle Kerne gleich- 
-zeitig an der Wasseraufnahme, so dass man sich fiiglich fragen kann, ob das von 
_JuncE angenommene verschiedenartige Verhalten der verschiedenen Kerne bei 
diesen Messungen zu Recht besteht. 
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Macue [10] bewies mit einer von ihm ersonnenen Vorrichtung, dass die durch 
adiabatische Expansion gebildeten Nebeltrépfchen sich ausserordentlich rasch auf- 
bauen und bei Kompression ebenso rasch wieder verdampfen und so einem Druck- 
wechsel bis zu 331/,; Hz zu folgen vermégen. Diese Beobachtung fiihrte ihn dazu, 
eine so schnelle Folge von Kondensation und Verdampfung anzunehmen, dass nach 
seiner Auffassung die Luftprobe sich nur im Zustande der Feuchtesattigung be- 
findet und somit keine Ubersadttigung in Erscheinung tritt. 

Da MacueE nur Nebel im Moment der gréssten Kompression bzw. der starksten 
Expansion beobachtete, geniigen seine Versuche nicht, zu entscheiden, ob nur eine 
Auslese oder die Gesamtheit der Kerne sich an der Nebelbildung beteiligen. Ausser- 
dem erhebt sich die Frage, wie rasch Nebelbildung und Druckanderung aufeinan- 
der folgen. 


2. Die Messapparatur 


Die zur Erfiillung der angedeuteten Aufgabe entwickelte Messapparatur ge- 
stattete, sowohl die adiabatische Zustandsanderung der Luftprobe als auch die 
Triibung eines die Luftprobe durchstrahlenden Lichtbiindels infolge Nebelbildung 
als Funktion der Zeit zu ermitteln. Damit kennen wir auch die Triibung (das heisst 
die Tropfenbildung) als Funktion der Zustandsanderung der Luftprobe. 

Das Schema der Messapparatur ist in Figur 2 dargestellt. Als Druckmessgerat 
diente ein Philips-Druckdifferenzgeber PR 6000 F, welcher mit einem Quecksilber- 
manometer geeicht worden war. Dieser wandelte eine Druckdifferenz in eine elek- 
trische Spannung um, welche der Druckdifferenz innerhalb + 2,5% linear folgte; 
die Spannung wurde mit einem Philips- Kathodenstrahloszillographen GM 5666 ge- 
messen. Die tiefste Eigenfrequenz des Gebers lag bei 190 Hz, so dass statische und 
dynamische Messungen bis etwa 100 Hz ohne nennenswerte Verzerrungen mdglich 


sind. Die Figuren 3 bzw. 4 zeigen Aufnahmen des Mess- bzw. Registrierteils der 
Apparatur. 


Messung der Extinktion 


Ein paralleles Strahlenbiindel wurde zwischen zwei im Expansionsgefass be- 
findlichen Spiegeln mehrfach reflektiert, so dass eine Weglange des Lichtes zwischen 
60 und 200 cm eingestellt werden konnte. Als Lichtquelle diente eine Gliihlampe 
mit 25 W Volleistung, welche aus Batterien gespeist wurde. Nach Durchlaufen des 
Rezipienten traf der Lichtstrahl die Vakuumphotozelle (Philips 90 CV), welche 
iiber einen Verstarker eine der Anderung der Lichtintensitat proportionale Span-_ 
nung an einen Dumont-Kathodenstrahloszillographen 304 A gab. Durch geeignete 
aed ist iiberpriift worden, dass das System Photozelle-Verstarker in dem 

ereic 


ae oo ee eee 

=7, { 

: 

linear arbeitete, wobei J, die Intensitat des ungetriibten, J diejenige des durch die 
Wolke getriibten Lichtes bedeutet. Die Bezugsspannung am Verstarker konnte 
beliebig eingestellt und auf diese Weise die ganze Schirmhéhe des Oszillographen- 
bildes ausgenutzt werden. A 


Um die Gleichzeitigkeit von Druck- und Extinktionsmessung zu gewalirleisteil 
erfolgte das Ingangsetzen der Anlage durch eine auf dem Schlagbolzenprinzip ar-_ 
beitende Vorrichtung. Diese schaltete zunachst iiber ein Verzogerungsrelais einen’ 
Multivibrator ein, welcher einen positiven und nach 10 m ; 


s einen negativen ne 
nungsstoss erzeugte. Der erste Sp . ie : 


annungsstoss schaltete die einmalige Zeitablenkung — 


' 


Vol. X, 1959 Kurze Mitteilungen — Brief Reports - Communications bréves 513 


AuslaBventi! 


Gebergruppe 


. 
| 


Rezipient 


PR 6000F 


Registriergruppe 


Philips 6M 5666 Du Mont 304A 


Figur 2 
Prinzipschaltbild der Messapparatur. 


Figur 3 
= Aufnahme der Registriergruppe. 


a 

As 
y 

J 
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Figur 4 


Aufnahme der Gebergruppe. 


der beiden Kathodenstrahloszillographen ein, wahrend der zweite in deren Z-Ab- 
lenkung eingespeist wurde; dadurch wird nach 10 ms Laufzeit auf den Strahlwegen 
beider Oszillographen gleichzeitig eine Lichtmarke angebracht. Zur Festlegung des 


zeitlichen Verlaufs der Strahlablenkung ist beiden Oszillographen eine bekannte 


Frequenz iiberlagert und sowohl Strahlverlauf wie Eichfrequenz fotografiert wor- 
den. Zur Erzeugung der Eichfrequenz diente ein Philips-RC-Generator GM 2315; 
in einigen Fallen verwendeten wir dazu die Eichspannung der Kathodenstrahl- 
oszillographen. 

Weniger als 0,1 ms nach dem Einschalten des Multivibrators 6ffnete der Schlag- 
bolzen das Offnungsventil. Am Kathodenstrahloszillographen konnte kontrolliert 
werden, ob der Beginn der Strahlablenkung genau mit dem Beginn der Druck- 
anderung zusammenfiel. Die Innehaltung dieser Bedingung ist fiir die Auswertung 
der Aufnahmen wichtig, da die Zeitablenkung am Dumont-Kathodenstrahloszillo- 
graphen (fiir das Schreiben der Extinktion) in den ersten 10 ms nicht linear ist. 


Das Expansionsgefass 


Als Expansionsgefass wurde der in Figur 5 skizzierte Zylinder verwendet, mit 
welchem das Manometer, das Offnungsventil (tiber eine kurze Schlauchleitung zur 
Vermeidung der Ubertragung von Erschiitterungen beim Offnen), die Pumpe (iiber 
einen Absperrhahn) und der Druckgeber verbunden waren. Weiter befand sich im 
Zentrum des Expansionsgefasses ein Thermoelement (Létstelle leichter als 10-3 g), 
dessen thermische Tragheit, verglichen mit der Einstellzeit des Siemens-Licht- 
markengalvanometers, vernachlassigbar klein war. Im Deckel waren zwei Glaser 
fur den Ein- und Austritt des Lichtstrahles angebracht. Komplanar mit Deckel und 
Grundflache sind zwei Spiegel fiir die Mehrfachreflexion des Lichtstrahles zur Ver- 
grosserung der Weglange des Lichtes befestigt. Bei Untersuchungen der Wasser- 
dampfkondensation mussten die Spiegel und Glaser gegen Beschlagen durch einen 
Anstrich von «Permaklar» geschiitzt werden. Dieser Schutz geniigte, um die Inten- 
sitatsanderung des Lichtes infolge Wasserabscheidung an besagten Oberflachen 
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< 1%/min zu halten, welcher Wert, verglichen mit der durch Nebelbildung hervor- 
gerufenen Intensitatsanderung des Strahles, die auch bei den eaaites Expan- 
sionen > 1°/,,/10-* min betragt, von vernachlassigbarer Grésse war. Bei Utes: 
suchungen mit Benzin und Azeton waren Schutzmassnahmen nicht notwendi 
weil ein Beschlagen der Glaser und Spiegel nicht beobachtet werden konnte Ba 
Auslassventil wurde so dimensioniert, dass es durch Einbau von Blenden méglich 
war, den Druckausgleich iiber eine Druckanderung von 190 mm Hg (stets bezogen 
auf p*) wahlbar innerhalb 300 ms bis 6 s vorzunehmen. 


30 cm 
Figur 5 


Skizze des Expansionsgefasses. 


3. Durchfiihrung der Experimente und die Messergebnisse 


Der zylindrische Teil des Rezipienten wurde (von Offnungen fiir die Verbindung 
abgesehen) mit Fliesspapier ausgekleidet. Dieses ist mit der zu kondensierenden 
Fliissigkeit benetzt worden. Auch der Boden des Rezipienten wurde mit der zu 
untersuchenden Fliissigkeit belegt. Hierauf ist im verschlossenen, normale Zimmer- 
luft enthaltenden Rezipienten mit befeuchteter Luft ein Uberdruck von 190 mm Hg 
hergestellt worden. Bei Untersuchungen mit Benzin und Azeton wurde die Kom- 
pressionsluft vor Eintritt in den Rezipienten iiber eine Kiihlfalle mit — 14°C 
geleitet, um stérende Einfliisse des Wasserdampfes bei der Kondensation auszu- 
schalten. Nach Erreichen der gewiinschten Kompression erfolgte jedesmal eine 
_ Priifung dariiber, ob die Anzeige des Druckgebers der Eichkurve entsprach. Danach 
wurde abgewartet, bis die Erwarmung infolge Kompression sich im Rezipienten bis 
auf mindestens 0,1°C ausgeglichen hatte, was einige Minuten bendtigte. 

Unter der Annahme, dass mit der gleichen Geschwindigkeit auch die Sattigung 
_ angendhert wurde, kénnen wir voraussetzen, dass unmittelbaz vor Beginn der Ex- 
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pansion eine relative Feuchte von etwa 99% geherrscht hat. Ein volliger Tempe- 
raturausgleich bendtigte mindestens eine Stunde. Eine so lange Wartezeit — 
sich aber ungiinstig, weil sich in diesem Falle bereits eine erhebliche ee oe 
Kerndichte im Rezipienten zeigte. Von wenigen Versuchen abgesehen, wurde aher 
der Nachteil eines geringen Defizits bis zum Erreichen der Sattigung in Kauf 
eikeenee und die Registrierung an den beiden Kathodenstrahloszillographen 
sind in der oben beschriebenen Weise ausgefiihrt worden. Eine grossere Expansions- 
zeit als 6s stellte sich als unzweckmassig heraus, da nach etwa 10 Ss die Warme- 
iibergainge die Messungen zu beeinflussen begannen. Fir das Verstandnis der Bil- 
dung natiirlicher Wolken waren aber gerade kleinere Expansionsgeschwindigkeiten 
em Interesse gewesen. 

ee 6 bis 9 ne eine Auslese aus den angefertigten Aufnahmen. Und 
zwar ist in Figur 6 der Verlauf von Extinktion und Druckabnahme bei einer 
schnellen Expansion, in Figur 7 bei einer langsamen Expansion von wasserdampf- 
haltiger Luft wiedergegeben. In Figur 8 ist eine langsame Expansion von benzin- 
haltiger, in Figur 9 eine schnelle Expansion von azetonhaltiger Luft gezeigt. Die 
verwendete Eichfrequenz war bei allen vier Aufnahmen 50 Hz; 190: Lore Druck- 
anderung ergaben auf dem Oszillographen eine Anderung von 46 mm; bei den 
Figuren 6, 7 und 9 entsprach die Intensitatsdnderung 0 < J < J, 15cm, bei der 
Figur 8 25 cm Strahlablenkung. 


Figur 6 
Schnelle Expansion wasserdampfhaltiger Luft: 
a zeitlicher Verlauf der Druckanderung; 6 zeitlicher Verlauf der Extinktion. 


+75 + ‘Ae 25 


Figur 7 


=" Langsame Expansion wasserdampfhaltiger Luft: 
a zeitlicher Verlauf der Druckinderung; 6 zeitlicher Verlauf der Extinktion. 


ig Gg mm ML ae 
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Figur 8 
Langsame Expansion benzindampfhaltiger Luft: 
a zeitlicher Verlauf der Druckanderung; b zeitlicher Verlauf der Extinktion. 


175 25 
tem 
125 25; 


a b 
Figur 9 


ra 


Schnelle Expansion azetondampfhaltiger Luft: 
a zeitlicher Verlauf der Druckanderung; 6 zeitlicher Verlauf der Extinktion. 


Figur 6 demonstriert besonders deutlich, wie schnell die Kondensation erfolgt. 
Nach 10 ms ist die Lichtintensitat an der Photozelle bereits um 10%, nach 20 ms 
um 33% reduziert. Die Tropfenbildung folgt der Druckaénderung so rasch, dass 
eine zeitliche Verschiebung beider Gréssen nicht mehr feststellbar ist. Die Schreib- 
breite des Strahles entspricht dem zeitlichen Auflésungsvermégen der Aufnahme 
und liegt bei etwa 1 ms. Dies ist genau die gleiche Zeit, die beim Ablauf der Druck- 
anderung ben6tigt worden ist, um die Feuchtigkeit im Innern des Gefadsses von 
99% auf Sdattigungsfeuchte zu steigern. Nach 10ms war das Verhaltnis 
dp/p* = 0,0165 (mit p* = 920 Torr), nach 20 ms war dp/p* = 0,03. 

Die Druckanderung, mit welcher der in Figur 7 gezeigte Vorgang verlief, war 
grossenordnungsmassig 100mal langsamer als bei der Figur 6. Hier erkennen wir 
deutlich, dass innerhalb der ersten 50 ms die Zunahme der Kern- bzw. Trépfchen- 
grdsse so gering blieb, dass sie keine Triibung des Lichtes bewirkte, wobei das Ver- 
haltnis dp/p* in dieser Zeit etwa 0,005 betrug. Unter der Voraussetzung, dass die 
Zustandsanderung bis zum Erreichen der Sattigungsfeuchte trockenadiabatisch 
verlauft, berechnet man, dass die Druckanderung gerade ausreichend war, um die 
relative Feuchte von 99,4% auf Sattigungsfeuchte zu steigern. Nach weiteren 
120 ms ist dp/p* ~ 0,008 und die Lichtintensitaét an der Photozelle um etwa 22% 
geschwacht. Ein Vergleich mit Figur 6 lasst erkennen, dass bei der langsamen 
Expansion die Triibung etwa viermal so rasch erfolgt wie bei der schnellen Ex- 


pansion. : 


ZAMP X/33 
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Aus diesen Beobachtungen an Wassertropfen kann man zweierlei schliessen: 

1. Dass es lediglich apparativ bedingt war, wenn bei den Versuchen von MAcHE 
die oberste Frequenz, bei welcher man eine Nebelbildung noch beobachten konnte, 
bei 33!/, Hz lag. Aus der oben gefundenen Zeitdifferenz, nach welcher bei schnellsten 
Expansionen bereits eine deutliche Triibung festzustellen ist, ergibt sich, dass die 
Nebelbildung Druckanderungen bis etwa 200 Hz zu folgen vermocht hatte, wenn 
anstelle einer einmaligen Expansion ein mit gleicher Geschwindigkeit periodisch 
veranderlicher Druck angewendet worden ware. 

2. Da bei den Figuren 6 und 7 etwa die gleiche Anzahl an Kernen mit einem 
kleinen Scholzschen Kondensationskernzahler (8000—9000/cm’) gezahlt worden ist, 
kann man aus der verhdltnismdssig langsameren Ausbildung der Extinktion mit 
zunehmender Expansionsgeschwindigkeit ersehen, dass bei sehr grossen Expan- 
sionsgeschwindigkeiten der Ablauf der Kondensation durch Diffusionsvorgange 
beherrscht und begrenzt wird. 

Eine Méglichkeit, die Aufnahmen quantitativ auszuwerten, wird an anderer 
Stelle beschrieben werden. 

Bei den Versuchen mit Azeton und Benzin ist im wesentlichen das gleiche Ver- 
halten wie bei den Versuchen mit Wasser beobachtet worden; die durch Benzin- 
wolken bewirkte Extinktion scheint etwas schneller, jene von Azetonnebeln hervor- 
gerufene Triibung etwas weniger rasch zu erfolgen als die Triibung bei Wasserwolken. 

Die oben beschriebenen Befunde sind im wesentlichen bei allen Aufnahmen 
festgestellt worden. Damit erkennen wir als wichtigstes Ergebnis der vorgenom- 
menen Messungen: 

Bei langsamen Expansionen hinkt die Wolkenbildung hinter der Druckanderung 


um eine kleine Zeitdifferenz von etwa 50 ms nach. Diese selbe Zeit wird jedoch ~ 


grossenordnungsmassig gerade beansprucht, um die Feuchte vom Ausgangswert 
auf den Sattigungswert zu bringen. Bei schnellen Expansionen wird dagegen der 
Sattigungswert so schnell erreicht, dass im wesentlichen keine Aussagen tiber den 
Ablauf der Kondensation bis zur Sattigung mehr méglich sind, da das zeitliche 
Auflosungsvermégen des benutzten Registrierverfahrens dazu nicht geniigte. Auch 
beeintrachtigte bei schnellen Expansionen die Tragheit des Diffusionsvorganges 
bereits die Kondensationsgeschwindigkeit. 

Der Firma Philips AG, Ziirich, méchten wir fiir die kostenfreie Uberlassung eines 
Druckdifferenzgebers, eines RC-Generators sowie eines Kathodenstrahloszillogra- 
phen zur Durchfiihrung der Versuche danken. 

Fiir das fordernde Interesse, welches Herr Professor Dr. R. SANGER, ETH, 
Ziirich, und der Leiter des Osservatorio Ticinese, Herr J.C. THams, dieser Arbeit 


entgegengebracht haben, erlauben wir uns, an dieser Stelle unseren besonderen 
Dank auszusprechen. 
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Summary 


In accordance with KéuHLER's theory, considerable supersaturation is needed to 
get the entire body of condensation nuclei to grow into droplets; the degree of 
saturation is determined by the size of the individual nuclei and their solubility in 
the condensed liquid. Since it is known that in expansion meters all the atmos- 
pheric particles which count as condensation nuclei do in fact grow into droplets, 
it may be assumed that the designated supersaturation in the expansion chambers 
is achieved for at least a short space of time. The formation of the theoretically 
correct degree of supersaturation requires, however, that once the moisture satura- 
tion has been exceeded condensation should markedly lag behind pressure variation, 
in order to allow time for the requisite supersaturation to develop. 

Experiments designed to test this question have established that condensation 
follows extraordinarily quickly on variation of pressure after moisture saturation 
has been exceeded. Condensation is also virtually uninfluenced by the choice of 
vapour used. At very high expansion speeds the condensation process is governed 
by the operative laws of diffusion. Measurements indicate that the degrees of super- 
saturation occurring in the expansion chambers lie considerably below the level of 
critical saturation even for the smallest nuclei. 


(Eingegangen: 17. Februar 1959.) 


Untersuchungen iiber die partielle Oxydation von pulverformigem 
Kupferoxydul zu Kupferoxyd 


Von PETER BRUESCH, Ziirich?) 


1. Problemstellung 


Die nachstehenden Untersuchungen iiber die partielle Oxydation von pulver- 
férmigen Cu,O zu CuO sind im Zusammenhang mit Messungen iiber die Eisbil- 
dungsaktivitat durchgefiihrt worden. Sowohl Cu,O als auch CuO sind in der Lite- 
ratur®) als Eisbildungskerne bekannt, das heisst, sie vermégen bei einer verhaltnis- 
missig hohen Temperatur die Vereisung der Wassertropfchen eines unterkthlten 
Nebels einzuleiten. Es wird vermutet, dass die Oberflache der Eisbildungskerne 
entscheidend an der Eiskeimwirksamkeit beteiligt ist. Um die Rolle der Oberflache 
abzuklaren, kann zum Beispiel versucht werden, die Oberflachen der Kerne auf 
méglichst iibersichtliche Weise zu stéren und den Einfluss dieser Stérung aut die 
Eiskeimaktivitat festzustellen. Die naheliegendste Stérung besteht in der partiellen 
Oxydation des Cu,O zu CuO. In der vorliegenden Arbeit wird das Ziel vertolgt, 
jedes einzelne Teilchen des pulverformigen Cu,O mit einer zusammenhangenden 
CuO-Schicht messbarer Dicke zu umgeben. 


1) Forschungsstelle fiir Physik der Atmosphare an der ETH. , 
2) R. SANGER, The Mechanism of Ice-Forming Nucleability, Bull. Obs. Puy de Dome 76 (1958). 
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2. Die partielle Oxydation 


Voraussetzung fiir die homogene partielle Oxydation ist die Dispergierung des 
Cu,O-Pulvers in die das Pulver zusammensetzenden Teilchen. Vorversuche haben 
gezeigt, dass gréssere Konglomerate nicht homogen oxydiert werden kénnen, da 
sich der Sauerstoff nur an der Oxydation der dusseren Teilchen beteiligt, wahrend 
die inneren Teilchen des Konglomerates im Oxydationsschatten liegen. Die Disper- 
gierung der Teilchen erfolgt mit einer sogenannten Windsichtungsapparatur: Ein 
senkrecht aufsteigender Luftstrom fiihrt die Teilchen durch ein konisches Glas- 
gefass, so dass die Luftgeschwindigkeit von unten nach oben stetig abnimmt. Die 
Luftgeschwindigkeit ist so gewahlt worden, dass die die Sichtungsapparatur ver- 


lassenden Teilchen Durchmesser von 1 bis 2 u aufweisen. Der die Sichtungsappara- — 


tur verlassende Luftstrom fiihrt die dispergierten Teilchen durch einen zylindrischen 
Oxydationsofen von 100 cm Lange und 7 cm Durchmesser. Die partiell oxydierten 
Teilchen werden schliesslich auf einem Glasfilter aufgefangen. Alle Versuche sind 
mit konstanter Sauerstoffkonzentration (Sauerstoffgehalt der Luft) und konstanter, 
dusserst geringer Teilchenkonzentration durchgefiihrt. Die Luftgeschwindigkeit im 
Ofen betragt etwa 15 cm/min, so dass die Teilchen geniigend Zeit haben, die Maxi- 
maltemperatur im Ofen anzunehmen. (Die im folgenden angegebenen Temperaturen 
bedeuten die Maximaltemperaturen im Oxydationsofen.) Der fiir die partielle 
Oxydation interessante Temperaturbereich liegt fiir diese. Versuchsanordnung 
zwischen 330°C und 470°C. 


3. Die Analyse der partiell oxydierten Proben 


Bei der partiellen Oxydation von Cu,O zu CuO sind folgende 4 Falle denkbar: 
a) Jedes Cu,O-Teilchen ist von einer zusammenhangenden CuO-Schicht umgeben 

und daher vollstandig von dieser eingekapselt; 

b) die CuO-Schicht ist nicht zusammenhangend; das heisst, nur ein Teil der Ober- 
flache ist durch CuO zugedeckt; 

c) es entsteht ein Gemisch von reinem Cu,O- und reinen CuO-Teilchen; 

d) Kombinationen von a), b) und c). 

Nimmt man an, dass bei der partiellen Oxydation eine zusammenhangende 
CuO-Schicht entsteht (Fall a), deren Schichtdicke iiberall mindestens so gross ist 
wie die Eindringtiefe der Elektronenstrahlen, so zeigen die Elektronenbeugungs- 
diagramme dieser Praéparate nur die CuO-Linien, wahrend die Réntgendiagramme 
bei geniigend dicker CuO-Schicht sowohl CuO- als auch Cu,O-Linien aufweisen. 
Bei den Moglichkeiten b), c) und d) enthalten dagegen die Elektronenbeugungs- 
diagramme CuO- und Cu,O-Linien. Erst die Kombination von Réntgen- und Elek- 
tronenbeugungsuntersuchungen lassen also eindeutige Resultate erwarten. 

Um Aufschluss tiber das bei der partiellen Oxydation auftretende Mengenver- 
haltnis von Cu,O zu CuO zu erhalten, werden Mischungsproben hergestellt, begin- 
nend mit 1,5% Cu,O, 98,5% CuO bis 98,5% Cu,O, 1,5% CuO, insgesamt 15 Proben. 
Es zeigt sich, dass die Nachweisempfindlichkeit des Cu,O im CuO etwa 1,5%, die- 
jenige des CuO im Cu,O etwa 3% betragt. Aus den Intensitiitsverhaltnissen der 
Cu,0- und CuO-Linien, die aus den Diffraktometerdiagrammen bestimmt werden 
k6nnen, ergibt sich eine Eichkurve, die gestattet, zu jedem Intensitatsverhdltnis 
das zugehérige Mischungsverhdltnis zu entnehmen. Auf diese Weise lasst sich aus 
den beobachteten Intensitaétsverhaltnissen der partiell oxydierten Proben deren 
Mengenverhaltnis zwischen CuO und Cu,O bestimmen. 

Rontgen- und Elektronenbeugungsuntersuchungen erfolgen direkt auf den auf 
den Glasfilter angelagerten partiell oxydierten Cu,O-Teilchen. Die gleichen Glas- 


Va. soe Ne AL 


a 
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filter konnen sowohl im Réntgendiffraktographén als auch im Elektronendiffrakto- 
graphen verwendet werden, was die Vergleichsmé@glichkeit wesentlich verbessert 

Die folgende Serie der Diffraktometerdiagramme (Figur 1) zeigt die Intensitats- 
verhaltnisse der partiell oxydierten Cu,O-Proben im Temperaturbereich zwischen 
320°C und 470°C. Fiir die Intensitatsvergleiche ziehen wir die beiden stairksten 
Cu,O-Linien (111) und (200) sowie die beiden starksten CuO-Linien (111) und (111) 
miteinander in Vergleich. Um genauere Intensitatsvergleiche zu gewinnen, wird der 
Winkelbereich zwischen 34° und 44° in oszillierender Bewegung zweimal abgetastet. 

Das Oxydationsdiagramm Nr. 1 der Figur 1 mit einer Oxydationstemperatur 
von 330°C zeigt nur die beiden Cu,O-Linien (111) und (200). Der prozentuale Anteil 
des CuO ist deshalb gemiass der Eichkurve kleiner als 3%. Das Oxydationsdiagramm 
Nr. 2 der Figur 1 mit einer Oxydationstemperatur von 350°C enthalt bereits Andeu- 
tungen der beiden staérksten CuO-Linien (111) und (111). Das zugehérige Intensi- 
tatsverhaltnis ergibt dabei einen CuO-Gehalt von etwa 3%. Wird die Oxydations- 
temperatur weiter gesteigert, so nimmt der CuO-Gehalt rasch zu. Bei 470°C sind 
keine Cu,O-Linien mehr nachweisbar. 

Die Entscheidung, ob es sich um den erwiinschten Fall a) oder um einen der 
3 andern Falle handelt, bringen nun die Elektronenbeugungsdiagramme. Das Stu- 
dium der Anlagerung des CuO auf den Cu,O-Teilchen ist vor allem fiir kleine pro- 
zentuale CuO-Gehalte interessant, das heisst im Temperaturbereich zwischen 320°C 
und 360°C. Die 6 Aufnahmen der Figur 2 zeigen die Elektronenbeugungs- 
diagramme im kritischen Temperaturbereich. 

Das erste Elektronenbeugungsdiagramm mit einer Oxydationstemperatur von 
320°C stellt das Diagramm des reinen Cu,O mit den folgenden Interferenzen dar: 


hkl: (110) (4111) (200) (211) (220) (311) 
(ey foot 100 37 1 Ay 17 
Dabei bedeuten /k/ die Millerschen Indizes der Netzebene hAk/, und I,,, relative 
Intensitaten. 

Das letzte der Elektronenbeugungsdiagramme mit einer Oxydationstemperatur 
von 360°C zeigt das Diagramm des reinen CuO mit den folgenden Interferenzen: 


hkl: (110) (111) (111) (200) (202) (113) 


| ee 100 96 30 25 20 
hkl: (022Y (311) (310) (220) (113) (311) (312) 
ees to 15 19 7 


Da die Linien (022), (311), (310), (220), (113) sehr nahe beieinander liegen, erscheinen 
sie als eine einzige diffuse Linie. In den Diagrammen Nr. 2 bis 5 werden mit stei- 
genden Temperaturen die Cu,O-Linien schwacher, die CuO-Linien starker. So 
macht sich im 2. Diagramm mit ¢ = 330°C bereits die CuO-Linie (111) bemerkbar 
sowie die aus den obengenannten 5 Komponenten bestehende Interferenzlinie. Das 
3. Diagramm mit ¢ = 338°C zeigt schon ziemlich deutlich die beiden starksten CuO- 
Linien (111) und (111), die aber noch nicht voneinander getrennt erkennbar sind, 
weil die zwischen ihnen liegende Cu,O-Linie (111) die Trennscharfe beeintrachtigt. 
Beim 5. Diagramm treten bereits alle CuO-Linien in Erscheinung, und als einzige 
Cu,O-Linie kann nur noch die Linie (200) erkannt werden. Ist die Oxydationstem- 
peratur hoher als 350°C (Diagramm Nr. 6 mit ¢ = 360°C), so weist das Beugungs- 
bild nur noch die CuO-Linien auf, wahrend das entsprechende Rontgendiagramm 
der gleichen Probe (Diagramm Nr. 3, Figur 1) immer noch CuO- und Cu,O-Linien 
enthalt. Aus der Eichkurve entnehmen wir, dass die Oxydationstemperatur 


ZAMP 
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470°C 


Diagramm Nr. 9, ¢ 


= 460°C 
x0) 


jgur 1 
igramme der partiell oxydierten Proben. 


98% Cu0/2% Cu 


Diagramm Nr. 8, ¢ 


d 


Rontgendi 


= 44°C 
u,0 


7; % 
% C 


0% CuO/20 


8 


Diagramm Nr. 
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bug 
es aN (22018 
220 022 311,310 
271 200 111 


200 


71 
710 


177 


Diagramm Nr. 1, ¢ = 320°C Diagramm Nr. 2, ¢ = 330°C 


Cuz0 


317 

31 220,113 ee 

220 022,811,310 

200. m1 

m 
Diagramm Nr. 4, t = 345°C 
Cug0 Cu 0 Cu 0 
{371 . 
(312 ! 
220 | (377,312 
113 (220,113 
022 |022311.310 
“ie 713 
370 : 
113 ate 
200 202 200 
200 717 
171 117 
711 110 
Diagramm Nr. 5, ¢ = 350°C Diagramm Nr. 6, 1= 360°C 
Figur 2 


Elektronenbeugungsdiagramme der partiell oxydierten Proben im kritischen Temperaturbereich. 
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t — 360°C zu einem CuO-Gehalt von rund 7% fiihrt. Wir stellen damit fest, dass der 
gefundene CuO-Gehalt von 7% aus einer zusammenhangenden Schicht bestehen 
muss, welche jedes Cu,O-Partikel umschliesst. Pes 

Ist der CuO-Gehalt hingegen kleiner als 7%, so kann nicht mit Sicherheit auf 
eine zusammenhangende CuO-Schicht geschlossen werden, denn in diesem Fall 
zeigen die Elektronenbeugungsdiagramme sowohl CuO- als auch Cu,O-Linien, Dabei 
konnen einerseits die Cu,O-Interferenzen von freiliegenden Cu,O-Partien an der 
Oberflache herriihren oder andererseits Cu,O-Linien darstellen, deren zugehorige 
Schichtdicke kleiner ist als die Eindringtiefe der Elektronenstrahlen. In diesem 
Fall kann der darunterliegende Cu,O-Kern ebenfalls zu Cu,O-Interferenzen Anlass 
geben. (Der Begriff einer einheitlichen Schichtdicke ist mit Vorsicht anzuwenden, 
denn es besteht die Méglichkeit, dass bestimmte Netzebenen fiir den Aufbau der 
CuO-Schichten besonders giinstig sind, andere weniger giinstig. Es kann sich also 
nur um eine mittlere Schichtdicke handeln.) 

Eine approximative Berechnung zeigt, dass die mittlere Schichtdicke exponen-_ 
tiell mit der Oxydationstemperatur wachst. Die diinnste, noch nachweisbar zu-~ 
sammenhangende Oxydationsschicht ergibt auf Grund einer approximativen Be- 
rechnung eine Schichtdicke von etwa 120 A. 3 


Es ist mir eine angenehme Pflicht, Herrn Prof. Dr. R. SANGER fiir die Uber- 
tragung der Aufgabe und das grosse Interesse am Fortgang der Arbeit an diesem 
Stelle herzlich zu danken. Ebenfalls Dank schulde ich Herrn U. Katz, Mitarbeiter 
am Laboratorium fiir Atmospharenphysik der ETH, ftir seine wertvollen Rat- 
schlage, Herrn Prof. Dr. W. Epprecut, Abteilung fiir industrielle Forschung, fiir” 
seine Anregungen und Uberlassung der Elektronenbeugungsapparatur, Herrn— 
Dr. S. STEINEMANN, Laboratoire Suisse de Recherches Horlogéres, fiir die erstem™ 
Elektronenbeugungsaufnahmen, sowie Herrn Dr. B. BOHLEN, Technisch-chemi- 
sches Laboratorium der ETH fiir die ersten Versuche zur partiellen Oxydatio: 


und Herrn Dr. W. Jurzr, Gaslaboratorium der EMPA, fiir die Uberlassung der 
Windsichtungsapparatur. 


Summary 


X-ray and electron-diffraction tests show that, with the aid of a combined 
oxidizing oven and air sifting apparatus, it is possible to achieve partial oxidation 
of the Cu,O particles in a powder in such a way that each Cu,O particle is encased 
within a composite CuO layer. The thinnest CuO layer which can still be shown to be 
composite has a mean thickness of about 120 A according to an approximate 
method of calculation. The procedures adopted for this experiment cannot deter- 
mine whether still thinner layers are also composite. The thickness of the CuO 
layers increases exponentially with the oxidation temperature. This work has 


been undertaken in conjunction with research into the ice-forming nucleability 
of substances. 


(Eingegangen: 8. Dezember 1958.) 
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Ermittlung der Wurzeln einer Gleichung 4. Grades durch Aufspalten 
in quadratische Faktoren 


Von SILvio CHRISTELLER, Basel 


Eine Gleichung 4. Grades mit reellen Koeffizienten der Gestalt 
Mi+tas+be+c=0 
lasst sich immer in zwei reelle quadratische Faktoren zerlegen: - 
M+axvtitbe+ce= (+544 p) (2 —sxrtqg=0. (1) 


Diese Aufspaltung wird im folgenden durchgefiihrt nach einem Probierverfahren, 
das eine Spezialisierung eines von ZURMUHL [1]!) und BuckincHam [2] benutzten 
Verfahrens ist, wobei aber zur Verbesserung der Naherungswerte die Methode von 
NeEwTon-Rapuson benutzt wird. 

Aus (1) folgt durch Koeffizientenvergleich unmittelbar 


ptq—st=a, (2) 
S=?) =: (3) 
Page. (4) 
Man wahlt nun f und berechnet aus (4) g = c/p und dann 
= b a b 
: Pgh mT ae 


aus (3). Schliesslich kontrolliert man, ob die Relation (2) erfiillt ist, das heisst, man 
prift, ob die Testfunktion 


bo 


Aes 3 | b 
eet lees) (5) 


gleich 0 ist?). 
Differenziert man T(p) nach #, so ergibt sich 
é b 22 cop+p 
T'(p) =1 2 ( ) 6 
= pe \cp—p) pb clp— 2 
‘Die Berechnung der Ableitung macht hier gar keine Schwierigkeiten und ergibt 
sich ohne wesentliche Mehrarbeit*). 


1) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 527. 

2) Auch HARTREE geht vom Ansatz (1) bis (4) aus, wahlt aber s und verbessert durch inverse 
Interpolation; vgl. [3], S. 220-221. 

3) Das in [1] angegebene Verfahren geht von drei Werten einer Testfunktion fiir drei aqui- 
distante Werte der unabhangigen Variablen aus. Verbessert wird durch quadratisches Eingabeln, 
da die Berechnung der Ableitung der Testfunktion zu umstandlich ist. In [2] wird dagegen linear 


interpoliert. 


' ; 
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Durch grobes Probieren findet man immer einen Naherungswert fiir p in einem — 
der Intervalle (—k, 0) oder (0, #), wo zur Vereinfachung k = || c| gesetzt wurde. 
Ist namlich c > 0, dann folgt aus (5) 


T(+0)=+00, T(k) =—oco 


und aus (6) 
TY (py 10 (Op RY 


In (0, k) existiert daher genau ein py mit T(po) = 0. 
Im Fallec < 0 hat man zwei Falle zu unterscheiden. Zunachst beachte man, dass 


T(+0) =—o, T(-0)=+o0, T(—co) = —oo 
ist. Ferner gilt nach (4) und (5) 
Defoe (2) Sens 


a) Fiir T(k) <0 existiert in (—k, 0) genau ein fy mit T(po) = 0 wegen 
T’(p) > O fir —k < p<0. 

b) Fiir T(k) > 0 gibt es mindestens ein p, in (0, k) mit T(p 9) = 0. Unter Um-— 
standen ist T’(p,) sehr klein oder 0. Dann kann man aber p in (— oo, —k) wahlen. 

Die ganze Rechnung erfolgt am besten nach einem Schema, wie im folgenden 
Beispiel (vgl. [1], S. 61): 


a4 — 38+ 2024+ 24% —-7=0. (7)a 


Setzt man * = y + 0,75, so erhalt man ‘ 
“ 


yt — 1,375 y2 + 1,625 y — 5,32421875 = 0. (8) 


Somit ist a = — 1,375, b = 1,625 und c = — 5,32421875. Die beiden Kolonnen fiir” 
p = 2 und # = 1,75 wurden mit Hilfe eines Rechenschiebers gerechnet. i 


‘ 


| p | 2 | 1,75 | 1,76 1,762029 | 1,76203127 | 


— 2,662 | — 3,055 | — 3,025124 — 3,021 64081} — 3,02163692 | 
— 0,662 | — 1,305 | — 1,265124 — 1,25961181] — 1,25960565 | 
@-p — 4,662 | — 4,805 | — 4,785124 — 4,783 66981 | — 4,783 66819 | 

— 0,339594 — 0,339697 35 | — 0,339697 47 | 
0,115 324 0,115 39428 


ee 
me) 
4 
SS 


S. 
@ 0,121 0,114 


® 

@® 

® 

T=@-®-—a 0,682 | — 0,044 | — 0,005448 | — 0,00000609 
@ |@/p —1,718820 |-—1,71486440 
@ | 2/p 1,136 363 1,13505509 
@/@_ 0,264 386 0,263 314.95 
| © |©@® 0,034647 0,034.488 68 
/® | T’=1-@-@ 2,684173 2,68037572 
® | T/T’ = ©/® — 2,029 - 10-8] — 2,2721-10-8 


Nach zwei Schritten hat man p = p) = 1,76203127 auf 8 Dezimalen bestimmt. 
wobei der Fehler héchstens 10-8 betragt. | 


Damit ist (8) aufgespalten in 


(y? — 0,33969747 y + 1,762031 27) (y* + 0,33969747 y — 3,021 63692) = 0. 
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Fiir die Wurzeln findet man 


¥1= 0,16984873 + 1,316503967, 
Yo= 0,169848-73 — 1,31650396i, 
Ys = —1,916 41264, 
Y,=  1,57671518, 


und damit fiir die Wurzeln von (7) 


%,=  0,91984873 + 1316503967, 
%,= 0,91984873 — 1,31650396i, 
H, = —1,16641264, 
%_= 2,32671518. 


Die Summe der Wurzeln (7) betragt 3,000 00000 und ihr Produkt — 6,99999986. 
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Summary 


In order to determine the roots of an algebraic equation of the fourth degree 
with real coefficients and removed cubic term a method of splitting up in two 
quadratic factors is used. The method follows in principle a procedure used, e. g. by 
ZURMUHL and BUCKINGHAM to determine the coefficients of the quadratic factors, 
but permits the use of NEwron-Rapuson’s method instead of linear or quadratic 
interpolation. 


(Eingegangen: 13. Marz 1959.) 


Varia — Miscellaneous — Divers 


ILMAC - Internationale Fachmesse und Fachtagungen tiber 
Laboratoriumstechnik, Messtechnik und Automatik in der Chemie 
Basel, 10.-15. November 1959 


In Basel findet vom 10. bis 15. November 1959 die ILMAC, internationale 

_Fachmesse und Fachtagungen iiber Laboratoriumstechnik, Messtechnik und Auto- 

-matik in der Chemie, statt. Die Fachtagungen werden vom Schweizerischen Che- 

mikerverband und der Schweizerischen Gesellschaft fiir Automatik veranstaltet 

und bilden zugleich die 22. Veranstaltung der Europaischen Féderation fiir Chemie- 

_Ingenieur-Wesen. Die Fachmesse wird von den Organen der Schweizer Muster- 
messe durchgefiihrt. 
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Im Programm der Fachtagungen sind insgesamt 30 Hauptreferate vorgesehen. 
Es konnten dafiir namhafte Vertreter der Wissenschaft und der Praxis aus Belgien, 
Deutschland, England, Frankreich, Holland, den USA und der Schweiz gewonnen 
werden. Erganzend zum eigentlichen Tagungsprogramm sind noch Referate und 
Demonstrationen der an der Fachmesse beteiligten Firmen vorgesehen. 

Die Tagung fiir Laboratoriums- und Messtechnik (10. und 11. November 1959) 
ist auf das Gebiet der Chemie begrenzt mit Hauptreferaten iiber elektrische Metho- 
den, optische Methoden, Trennverfahren und Isotopenmethoden. Das Programm 
der anschliessenden Tagung iiber Automatik umfasst die technischen Fragen der 
Automation in der Verfahrenstechnik der chemischen Industrie und der Lebens- 
mittelindustrie, ferner die wirtschaftlichen Aspekte der Automation sowie die 
Datenverarbeitung in der Wirtschaft. 

Fiir weitere Auskunft wende man sich an das Sekretariat der ILMAC, Basel 21. 


errr rnin nner 
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Principles of Modern Physics. Von Rosert B. LEIGHTON (McGraw-Hill, — 
London 1959). 795 S., 207 Fig.; 97 sh. 

Das vorliegende Werk enthalt eine klare Einfiihrung in viele der modernsten 
und interessantesten Teilgebiete der neueren Physik. Als solche treten auf: die — 
Kernphysik mit den verschiedenen Kernmodellen, die Physik der Elementar- 
teilchen einschliesslich der Srangeness-Theorie von GELL-MANN und NISHIJIMA, 
das Problem der Paritatsverletzung bei schwachen Wechselwirkungen, die kos- 
mische Strahlung, die Energieerzeugung in Sternen und die Bildung der chemi- 
schen Elemente. 

Daneben sind auch die klassischen Theorien vertreten in Kapiteln tiber Atom- — 
und Molekiilspektren, die Bandertheorie des Festkérpers und die R6ntgenstrahlen. 

In Anbetracht der Mannigfaltigkeit des gebotenen Stoffes miissen die formalen 
mathematischen Entwicklungen auf ein Minimum beschrankt werden. Diese’ 
Beschrankung gelingt dem Verfasser oft in ausgezeichneter Weise. Die zwei ersten 
Kapitel sind der speziellen Relativitatstheorie und der Quantentheorie gewidmet. © 
Das erste dient wohl mehr der Repetition als der Einfiihrung. 

Musterhaft sind die knappen historischen Einfiihrungen und die Zusammen- 
fassungen am Ende einzelner Kapitel. Die Schreibweise ist immer klar, lebendig 
und anregend. Die vielen Ubungsaufgaben erlauben dem Leser stets eine Kon- 
trolle seines Verstandnisses. . 

Das Buch kann besonders auch dem Studenten aufs warmste empfohlen ~ 


werden. R. Jost 


Colloque National de Magnétisme (Commémoratif de l’ceuvre de PIERRE 
WEIss). (Centre National de la Recherche Scientifique, Paris 1958) 336 S. 
174 Abb.; ffrs. 4000.-. . 

Das Kolloquium zu Ehren von PIERRE WEtss, dem Vater des Ferromagnetis- : 
mus, vereinigte die franzésischen Magnetophysiker vom 8. bis 10. Juli 1957 in ~ 
Strassburg. Mit der leider unvermeidlichen einjahrigen Verzogerung ist nun der 
Tagungsbericht erschienen, der 26 Arbeiten wiedergibt. Ein Drittel des Buches 
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fiillen die drei Hauptvortrage, in welchen die befruchtende Wirkung der Weiss- 
schen Molekularfeldhypothese auf die weitere Entwicklung der Theorie des 
Magnetismus gewiirdigt wird. N&rx gibt in seinem Vortrag einen Uberblick iiber 
die Erfolge der Wei8schen Vorstellungen und ihre Vervollkommnung durch den 
Autor selber, der mit Hilfe des lokalen Molekularfeldes das Verhalten von ferri- 
magnetischen und antiferromagnetischen Substanzen erklaren konnte. Eine griind- 
liche Einfiihrung in das Gebiet der Wei8schen Bezirke vermittelt der Vortrag 
von GUILLAUD und VauTIER. Die theoretischen Erkenntnisse werden an iiber 
zwanzig Photographien von Bitterfiguren diskutiert. PERRIER zeigt in seiner 
Arbeit, welche mehr oder weniger eine Zusammenfassung seiner Forschungen 
auf diesem Gebiet darstellt, den Einfluss des Molekularfeldes auf die elektrischen 
Eigenschaften (Hall-Effekt, Magnetowiderstand) von Ferromagneten. Die iibri- 
gen 23 Arbeiten bringen zum Teil neue Messungen und Ergebnisse, zum Teil auch 
abschliessende Berichte iiber Untergebiete oder bestimmte Substanzgruppen. 
Einen breiten Raum nehmen die Untersuchungen an Verbindungen der seltenen 
Erden, einschliesslich der Gruppe der Granate, ein, wahrend.Messungen an den 
historischen Ferromagneten Eisen und Nickel seltener sind, da diese besonders 
prazis sein miissen, damit sie etwas Neues bieten kénnen und eine bessere Priifung 
existierender Theorien erlauben als die schon von WEIss bestimmten Kurven. 
Natiirlich fehlen auch nicht Arbeiten itiber Halogenide der Eisengruppe, iiber 
Ferrite und Resonanzeffekte. Der Bericht gewahrt somit einen umfassenden 
Einblick in das Forschungsprogramm der franzésischen Magnetiker, die zur 
Entwicklung dieses Zweiges der Physik schon immer einen wesentlichen Beitrag 
leisteten. F. HULLIGER 


Theoretical Mechanics. An Introduction to Mathematical Physics. 
Von J. S. Ames und F. D. MurnaGuan (Dover Publications, London and New 
York 1958). 462 S., 39 Fig.; $ 2.00. 

Das Werk beginnt mit einem ausfiihrlichen Kapitel tiber Vektoranalysis und 
behandelt dann der Reihe nach die Kinematik des starren Korpers, die Dynamik 
des Massenpunktes und die Kreiseltheorie. Schwingungen und Wellen werden 
relativ kurz, die analytische Dynamik ausfiihrlicher behandelt, und ein letztes 
Kapitel fiihrt in die Relativitatsmechanik ein. Das Buch kann als Einfiihrung in 
die héhere Mechanik bestens empfohlen werden. H. ZIEGLER 


Handbuch der Physik - Encyclopedia of Physics. Herausgegeben von 
S. Fricce, Bd. 6: Elastizitat und Plastizitét (Springer-Verlag, Berlin 1958). 642 S., 
254 Fig.; DM 145.-. 

Der Band umfasst in sechs Beitragen von angelsachsischen Autoren die wichtig- 
sten Kapitel der Kontinuumsmechanik fester Korper mit Ausnahme der nicht- 
linearen Elastizitatstheorie. 

Ein erster Artikel von I. N. SNEDDON und D. S. Berry, The Classical Theory 


of Elasticity, vermittelt ausgewahlte Ausschnitte aus diesem weitlaufigen Gebiet. 
Die Grundgleichungen werden fiir endliche Deformationen und in der fiir die Be- 


handlung in krummlinigen Koordinaten geeigneten Tensorschreibweise entwickelt ; 
der Rest des Beitrags beschrankt sich im wesentlichen auf die kartesische Darstel- 


‘lung der Elastizitatstheorie infinimaler Verzerrungen. Neben dem Torsionsproblem 
werden ebene Aufgaben ausfiihrlicher behandelt, vor allem mit der von russischen 
-Mathematikern entwickelten komplexen Methode. Kapitel iiber raumliche und 


dynamische Probleme sowie eine kurze Darstellung der thermoelastischen Grund- 
gleichungen beschliessen den Beitrag. 
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Ein Artikel von H. T. Jessop, Photoelasticity, gibt eine Ubersicht iiber Theorie — 
und Praxis der zwei- und dreidimensionalen spannungsoptischen Verfahren samt 
Anwendungsbeispielen. ; 

A.M. FreupENTHAL und Hitpa GEIRINGER teilen sich in den Artikel The 
Mathematical Theories of the Inelastic Continuum, FREUDENTHAL gibt eine allge- 
meine, auch thermodynamische Betrachtungen umfassende Einfihrung in die 
heute aktuellen Gebiete der Kontinuumsmechanik und stellt vor allem in physi- 
kalischer Sicht die grundlegenden Beziehungen der Viskoelastizitat, Plastizitat 
und Viskoplastizitat zusammen. Hirpa GEIRINGER schliesst eine ausgezeichnete 
Darstellung der Fragestellungen an, die sich um den Begriff des idealplastischen 
K6rpers gruppieren, und diskutiert ausfiihrlich die mathematischen Verfahren zur 
Behandlung vorwiegend ebener Probleme mit Anwendungen auf elastisch-ideal- 
plastische und starrplastische Korper. 

Ein weiterer Beitrag von M. REINER, Rheology, schliesst sich an die bekannten 
Veréffentlichungen des Autors an und orientiert in oft eigenwilliger Darstellung 
hauptsachlich iiber die Fragestellungen, die fiir das Grenzgebiet Fest—Flissig — 
typisch sind. 

Ein kurzer Artikel von G. R. Irwin, Fracture, schildert die Behandlung der 
Ausbreitung von Rissen im festen Kérper mit energetisch-statistischen Verfahren, — 
und schliesslich gibt A. M. FREUDENTHAL unter dem Titel Fatigue einen zusammen- 
fassenden Beitrag zum Problem der Ermiidung. H. ZIEGLER 


Dreidimensionale Uberschallprobleme der Gasdynamik. Von K.-R. — 
DoRFNER ([Ergebnisse der angewandten Mathematik, Heft 3] Springer-Verlag, 
Berlin 1957). 150 S., 44 Abb.; DM 26.-. ; 

Der Verfasser hat sich zum Ziel gesetzt, die Grundlagen der linearen Methoden — 
zur Berechnung dreidimensionaler stationarer Uberschallstromungen unter Ver- 
wendung der in vielen Einzelveréffentlichungen verstreuten Untersuchungen syste- 
matisch darzustellen. Es handelt sich dabei im wesentlichen um das allgemeine 
Singularitatenverfahren und um die Methode der Uberlagerung kegelsymmetrischer 
Strémungen. Beide Hauptprobleme werden behandelt, das heisst die Berechnung 
des Strémungsfeldes, insbesondere der Druckverteilung, bei vorgegebener Geo- 
metrie des Tragfliigels, und umgekehrt die Bestimmung der Geometrie bei vorge- 
gebener Druckverteilung auf dem Tragfliigel. : 

Dem Charakter der Sammlung entsprechend, ist die Darstellung des Stoffes — 
betont mathematisch gehalten. Das Buch enthalt weder numerische Beispiele oder 
Zahlentafeln noch Versuchsresultate, so dass seine praktische Verwendung nur dem _ 
mathematisch sehr gut geschulten Aerodynamiker zu empfehlen ist. : 

Interessant ist die konsequente Anwendung der Théorie des Distributions von if 
LAURENT SCHWARTZ zur Darstellung unstetiger Funktionen. P. DE HALLER 
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Principles of Geodynamics. Von A. E. ScHEIDEGGER (Springer-Verlag, . 
Berlin 1958). 280 S., 86 Fig.; DM 49.60. : 

Dieses Buch gibt einen Uberblick iiber die Geschichte, den Aufbau, die Eigen-_ 
schaften und die Form der Erde. Es werden vor allem die geologischen, ozeanogra- ' 
phischen und geophysikalischen Erfahrungstatsachen zusammengestellt und die | 
zu ihrer Erklarung aufgestellten Theorien besprochen. 

Das ist offensichtlich eine weitschichtige und umfangreiche Materie. SCHEID- . 
EGGER hat aber durch eine geschickte und mit der notwendigen Kritik gepaarten 


Auswahl die Grundziige der Geodynamik in einer Form dargestellt, die neu ist und 
Anerkennung verdient. ‘ 
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Das Buch ist in folgende Kapitel aufgeteilt: I. Phys wogvaphic and Geological Data 
Regarding the Earth. 11. Geophysical Data Regarding the Earth. 111. The Mechanics 
of Deformation. IV. Effects of the Rotation of the Earth. V. Continents and Oceans. 
VI. Orogenesis. VII. Dynamics of Faulting and Folding. VIII. Dynamics of Other 
Features. 

Von den gestreiften Theorien sind natiirlich einige reine Spekulationen, und es 
ist auffallend, wie im Grunde genommen wenig durch Tatsachen gesichert ist. 

Das Buch ist allen, die sich mit der Erde als Ganzes befassen, vor allem zur 
Orientierung, sehr zu empfehlen. Max WEBER 


Handbuch der Physik - Encyclopedia of Physics. Herausgegeben von 
S. FLUcGE, Band 51: Astrophysik IT: Sternaufbau (Springer-Verlag, Berlin 1958) 
B30'., 197 Fig.; DM 175.—. 

Dieser zweite astrophysikalische Band im Rahmen des Handbuches der Physik 
gibt eine umfassende Darstellung des Aufbaues der Sterne, wobei in den meisten 
Beitragen iiber den eigentlichen Gegenstand hinaus auch Probleme der Stern- 
atmospharen behandelt werden, jedoch mehr im Sinne einer Erganzung als einer 
Duplizierung von Band I (Sternoberflachen, Doppelsterne). 

Der Band wird durch eine knappe und tibersichtliche Darstellung der Prinzipien 
des inneren Aufbaues und der Energieerzeugung der Sterne eingeleitet (M. H. Wru- 
BEL). Die Beitrage iiber das Hertzsprung—Russell-Diagramm (H.C. Arp) und die 
Haufigkeit der Elemente im Kosmos (H. E. Sugess, H. C. UREy und L. H. ALrEr) 
stellen die wichtigsten beobachtungsmassigen Grundlagen fiir das Studium der 
Sternentwicklung zusammen, mit welchem sich der 150 Seiten umfassende Artikel 
von E. M. und G. BursipGe, der das Kernstiick des vorliegenden Bandes darstellt, 

_beschaftigt. Entstehung und Entwicklung der Sterne werden von den verschieden- 
artigsten Aspekten — Kontraktion, Massenabgabe und -aufnahme, Kernreaktionen, 
Entstehung und Umwandlung der Elemente, Hertzsprung—Russell-Diagramme — 

beleuchtet. Dass die Resultate dieser jungen Forschungsrichtung noch vielfach 
widersprechend sind, wundert uns nicht. Der Artikel gibt deshalb weniger als die 
iibrigen eine enzyklopddische Ubersicht, dafiir viele stimulierende Ausblicke und 
Moéglichkeiten fiir die kiinftige Forschung. Der Band schliesst mit einigen meist 
kurzen Beitragen iiber spezielle Sterne und Sternzustande wie die verdnderlichen 
Sterne, die Sternstabilitat, die magnetischen Sterne, weisse Zwerge, Novae und 
Supernovae. Die Ausfiihrungen iiber Pulsation und Stabilitat der Sterne (P. LEDoux) ° 
sind, verglichen mit den andern Beitragen, unverhaltnism4ssig breit; man wird dies 
jedoch nicht bedauern, da es den Nichtspezialisten anderweitig schwer fallen wiirde, 
sich in dieser spréden Materie einen Uberblick zu verschaffen. M. WALDMEIER 


Wellenmechanik. Vorlesung von Prof. Dr. W. Pautt, ausgearbeitet von FRITZ 
Heriacu und Hernz E. Knorprer (Verlag des Vereins der Mathematiker und 
‘Physiker an der Eidgendssischen Technischen Hochschule, Ziirich 1959). 134 S., 
33 Abb.; sFr. 19.-. 

Inhaltlich besteht der Hauptunterschied zwischen dieser Kursvorlesung tiber 
‘unrelativistische Wellenmechanik und dem unrelativistischen Teil von PaUvLis 
‘beriihmtem Handbuchartikel darin, dass hier die prinzipiellen Fragen knapper, die 
mathematischen Methoden dagegen ausfiihrlicher behandelt werden. So wird unter 
-anderem das Beispiel des harmonischen Oszillators in der Ebene (in Polarkoordi- 
‘naten) und des Wasserstoffatoms (diskretes und kontinuierliches Spektrum, Lésung 
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in spharischen und in parabolischen Koordinaten) durchgerechnet. Die dabei auf-_ 
tretenden mathematischen Funktionen werden eingehend diskutiert. Besonderes 
Gewicht wird auf die konfluente hypergeometrische Funktion und ihre Spezialfalle 
gelegt. Weiter wird eine ganze Reihe von Naherungsverfahren behandelt: Neben 
der zeitunabhangigen und der zeitabhangigen Stérungsrechnung die Sattelpunkts- 
methode, die WKB-Methode und bei den Stossvorgangen die Bornsche Naherung 
und die Methode der partiellen Wellen. Dagegen wird das oft sehr niitzliche Varia- 
tionsprinzip nicht erwahnt. 

Der besondere Stil und die Originalitat von PauLis Vorlesung ist bei der Aus- 
arbeitung erhalten geblieben. Der Stoff ist sehr iiberlegen und elegant dargestellt, 
doch stellt der knappe Text an den Leser gewisse Anspriiche. Hin und wieder 
fehlt eine Erklarung, sei es, dass sie vergessen wurde (wie zum Beispiel die Definition 
der Operatoren #, p), sei es, dass Kenntnisse vorausgesetzt wurden, die man von 
einem jungen Studenten eigentlich nicht erwarten kann (eine derartige Lucke : 
klafft im Vollstandigkeitsbeweis fiir die Hermiteschen Polynome). Dafiir treten die 
wesentlichen Punkte, weil oft unbeschwert von mehr nebensdchlichen didaktischen ~ 
Details, besonders deutlich zutage. 

Einem Anfanger kann die Durcharbeitung dieser Vorlesung durchaus empfohlen 
werden, nétigenfalls unter Heranziehung eines der bekannten Lehrbiicher; sie 
wird ihm zu einer iiberdurchschnittlichen Beherrschung der Wellenmechanik ver- 
helfen. Fiir den Kenner der Wellenmechanik ist dieses Buch, dank der meister- 
haften Darstellung der mathematischen Funktionen, eine wertvolle Erganzung zum 
Handbuchartikel. Die Ubersichtlichkeit und ein Sachregister machen es auch als 
Nachschlagewerk geeignet. A. THELLUNG 


Das Mikroskop. Seine wissenschaftlichen Grundlagen und seine Anwendung. 
Von A. EnrtnGHAUS und L. TRApp (B. G. Teuber Verlagsgesellschaft, Stuttgart — 
1958). 144 S., 94 Abb.; DM 10.80. if 

Dem Bearbeiter dieser Neuauflage, Dr. L. TRapp, Wissenschaftlicher Mitar- 
beiter der Firma Carl Zeiss, Oberkochen, ist es nicht nur gelungen, dem in Mikro- 
skopikerkreisen bereits bekannten Ehringhaus ein neues Gesicht zu geben, sondern” 
sie auch auf den heutigen Stand von Wissenschaft und Technik zu bringen. Trotz 
des grossen Bedienungskomforts moderner Mikroskope sind gewisse Kenntnisse 
uber Funktionen und Leistung fiir die nutzbringende Anwendung der Instrumente — 
unentbehrlich. * 

In leichtfasslicher Weise behandelt das Buch zuerst die wichtigsten Begriffe 
und Grundlagen iiber Lupe, einfaches Mikroskop und die Abbildungsgeometrie des : 
zusammengesetzten Mikroskopes, einschliesslich die verschiedenen Abbildungsfeh- 
ler, um nach der wellenoptischen Betrachtung der Abbildung zur mehr praktischen - 
Seite iiberzugehen. } 

Neben der Beschreibung der Einrichtungen und Handhabung des Mikroskopes — 
und seiner optischen Systeme ist ein ausfiihrlicher Hauptabschnitt den verschiede- | 
nen Beleuchtungs- und Beobachtungsverfahren (einschliesslich Phasenkontrastver- 
fahren) gewidmet. Folgende Kapitel befassen sich mit Hilfsapparaten zum Mikro-— 
skop, der binokularen Beobachtung, Messungen .an mikroskopischen Praparaten 
und Leistungspriifung des Mikroskopes. Der Schlussabschnitt des Buches gibt noch > 
einen Uberblick uber die Zurichtung mikroskopischer Praparate. | 

So vermittelt das vorliegende Biichlein besonders auch dem interessierten Laien: 
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Gesammelte mathematische Abhandlungen 


von Lupwic ScniArct (1814-1895) 


In drei Banden, herausgegeben vom Steiner-Schlafli-Komitee 
der Schweizerischen Naturforschenden Gesellschaft 


Band I 


392 Seiten, mit einem Portrat von L. Schlifli 
In Ganzleinenband Fr./DM 54.— (1950) 


Band II 
381 Seiten. In Ganzleinenband Fr./DM 54.— (1953) 


Band III +: 
402 Seiten mit 20 Figuren. In Ganzleinenband Fr./DM 57.— (1956) 


BIRKHAUSER VERLAG - BASEL UND STUTTGART 


Neuerscheinung 


SCHIESSLEHRE 


von WALTER Stutz, Oberst der Artillerie 
Chef der Sektion fiir Schiessversuche 
des Eidgenéssischen Militardepartements in Thun 


(1959) 510 Seiten mit 358 Figuren und 75 Tabellen, gebunden Fr./DM 28.50 


Die Schiesslehre will weder ein Lehrbuch der Ballistik noch eine Einfiihrung in die 
Ballistik sein; es gibt vielmehr eine leichtverstandliche Orientierung iiber das weite 
Gebiet der Schiesslehre und vermittelt die notwendigen Grundlagen zum Verstand- 
nis der Waffen und deren zweckmissigem Einsatz, es vermittelt einen Uberblick. 


Aus dem Inhalt: 


Pulver und Sprengstoffe, Vorginge in der Waffe beim Schuss, Bewegung des Ge- 
schosses nach Verlassen der Miindung, Schiessaufgaben, Stérungen der Flugbahn, 
Richtlehre, GeschoBstreuung, Treffwahrscheinlichkeit, Schiessbehelfe, Geschosse, 
Ziinder, Schussladungen, Geschosswirkung, Bestimmen der Wirkungselemente, 

Wirkungsschiessen, Gesichtspunkte fiir den Einsatz von Wafien : 
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